4 


4 


#4 


ne 


Tome 21 


No 


Juizzer 1960 


LE JOURNAL DE PHYSIQUE 


ET 


LE RADIUM 


—— 


MÉCANIQUE STATISTIQUE QUANTIQUE. 
OPÉRATEURS DENSITÉS ET GRANDEURS THERMODYNAMIQUES. 


PATTES VON, 
GC. E. N. de Saclay. 


Résumé. — Calcul dans le cas de l’équilibre de l’opérateur densité simple et de l’opérateur 
densité double en fonction de la température et de l’activité. L’opérateur de diffusion permet 
dans certains cas de mettre à part les grandeurs thermodynamiques. En particulier, ces grandeurs 
ne figurent pas dans l’expression de l’opérateur densité double en fonction de l'opérateur densité 
simple : ce qui explique que cette expression puisse être valable même hors de l’équilibre. 


Abstract. — The simple and double density operators, as functions of temperature and activity 


are computed for a system in equilibrium. 
separating out the thermodynamical variables. 


The scattering operator allows, in some cases, 
In particular, these variables do not appear in 


the expression of the double density operator, as a function of the simple density operator ; which 
explains the fact that this expression is still valid out of equilibrium. 


I. Programme. — Les notations utilisées 1c1 
sont celles du mémoire intitulé « Une méthode 
d’étude des corrélations dans les fluides quantiques 
en équilibre » [1]. 

Le paragraphe 6 de ce mémoire était consacré 
à une tentative d’évaluation de l’opérateur den- 
sité simple y, et de l’opérateur densité double u,; 
grâce à un développement de l’exponentielle de 
Boltzmann-Gibbs suivant le thème imaginé par 
Ursell. Le résultat n’était pas très satisfaisant 
parce que les expressions obtenues ne paraissaient 
pas fournir, sans retouches artificielles, l’hermt- 
ticité indispensable. 

Lee et Yang [2] et Païs et Uhlenbeck [3] ont 
utilisé un développement meilleur à cet égard. Ils 
l’ont appliqué seulement à la « grande fonction de 
partition ». Nous nous proposons ici d'appliquer 
ce type de développement — avec des notations 
plus maniables — au problème des opérateurs 
densités. Cette entreprise nous paraît utile parce 
qu’il nous semble que la connaissance des opéra- 
teurs densités est plus fondamentale que celle de 
la grande fonction de partition. Ces opérateurs 
expriment mieux ce qu'est le «grain » du fluide, 
ils permettent notamment d'aborder directement 
les cas où le fluide n’est pas homogène, ils éli- 
minent bien des problèmes relatifs au volume du 
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récipient, enfin ils se prêtent à des manipulations 
que nous examinerons plus loin. 

La méthode d’Ursell n’est pas la seule : elle est 
en concurrence avec la méthode de perturbation 
qui a fait l’objet d'importants travaux. Je ren- 
verrai particulièrement aux mémoires de Bloch et 
De Dominicis [4]. Sans chercher à porter un juge- 
ment sur cette rivalité, disons tout au moins que 
la méthode d’Ursell mérite d’être explorée jusqu’au 
bout de ses possibilités en raison de sa simplicité. 

Quelle que soit la méthode employée, les opé- 
rateurs densités se trouvent finalement exprimés 
en fonction de la dynamique du système, de la 
température et de l’activité absolue. Il convient, 
pour essayer d’aller plus loin, de rapprocher les 
développements obtenus des développements ana- 
logues de la statistique classique. 

En statistique classique, les problèmes de répar 
tition dans l’espace sont complètement séparés 
des problèmes de répartition en quantité de mou- 
vement. La répartition de la quantité de mouve- 
ment est donnée par la loi de Maxwell. Ce point 
acquis, il ne reste plus qu’à calculer les densités, 
notamment la densité simple et la densité double. 
La méthode d’Ursell a fourni depuis longtemps dans 
ce cas des développements en fonction des forces, 
de la température et de l’activité — pour cette 
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dernière à l’aide d’une série entière. Malgré la 
séparation des quantités de mouvement dans le 
cas classique, il y a une analogie évidente entre 
les développements classiques et les développe- 
ments quantiques. 

Or les développements classiques demandent 
un certain remaniement si l’on désire qu’ils dé- 
crivent clairement la réalité physique. D’abord on 
résout l'équation qui exprime la densité simple en 
fonction du champ appliqué, des forces inter- 
moléculaires, de la température et de l’activité, 
en prenant comme inconnue le champ appliqué. 
On trouve finalement une équation où figure d’un 
côté une expression très simple du champ appli- 
qué, de la température et de l’activité et dont 
l’autre membre est une fonctionnelle de la densité 
— qui en général varie de point en point — et de 
la force intermoléculaire. Ce nouveau formalisme 
présente deux avantages : le premier — sauf 
peut-être si la densité est très élevée — est que 
l'équation est locale, en ce sens que, le champ 
appliqué ne figurant que par sa valeur en un point 
déterminé, les valeurs de la densité qui inter- 
viennent ne concernent que le voisinage de ce 
point (nous n’envisageons ici pour simplifier que 
le cas des forces à courte portée). Le second est 
que, l'expression du champ appliqué en fonction 
de la densité n’étant pas linéaire, pour un champ 
appliqué donné il n’y a pas qu’une répartition en 
densité : c’est la clé de la théorie cinétique des 
changements d’état. de l’équilibre des phases 
distinctes, et de l’existence des états métastables. 

De même, on parvient à exprimer la densité 
double en fonction de la densité simple, des forces 
intermoléculaires et de la température en élimi- 
nant le champ appliqué et l’activité. On obtient 
ainsi une représentation claire des corrélations 
statiques à l’intérieur du fluide. Cette représenta- 
tion est commode pour le calcul des moyennes 
qui font intervenir les molécules deux à deux : 
énergie potentielle, pression intermoléculaire, fluc- 
tuations en densité, etc... 

En vue de ce qui va suivre, il n’est pas sans 
intérêt d’évoquer une transformation de la repré- 
sentation examinée à l'alinéa précédent qui, 
pour la théorie classique, seule, n’est peut-être 
qu’une curiosité : au lieu d’exprimer la densité 
double en fonction des forces intermoléculaires, 
on recherche la relation inverse de manière à 
exprimer l’énergie potentielie intermoléculaire en 
fonction de la densité double. On obtient une rela- 
tion du type suivant : 


énergie intermoléculaire fonctionnelle de la densité 
température = simple et de la densité double. 


Iei, il y a une séparation complète entre la dyna- 
mique et la température imposée, qui sont d’un 
côté, et les grandeurs statistiques qui sont de l’au- 
tre. Une relation encore plus dépouillée s’ob- 
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tient en exprimant la densité triple en fonction 
de la densité simple et de la densité double : 
la dynamique — aussi bien les interactions que 
le champ appliqué — la température, l’activité 
ont disparu. Il s’agit donc d’une relation statis- 
tique assez fondamentale. La fonctionnelle obte- 
nue comporte une suite de termes de complica- 
tion croissante dont le premier, s’il est seul retenu, 
exprime l’approximation de Kirkwood. 

Il est tentant de rechercher des relations du 
même genre dans la théorie quantique. Un exemple 
très simple en a été obtenu, mémoire [1], formule 
(20), à propos de l’approximation du champ self 
consistant. Ce résultat a ‘été largement généralisé 
— pour un fluide homogène, pour des interactions 
fortes, par la méthode de perturbation — par 
Balian et De Dominicis [5]. Un phénomène 
dégagé par ces auteurs et surprenant au premier 
abord est une séparation plus immédiate encore 
qu’en théorie classique des variables thermody- 
namiques, activité et température : ceci tient au 
rôle sous-jacent de la répartition de Maxwell 
dans la théorie classique. La séparation classique 
entre espace et quantité de mouvement est d’une 
grande commodité, mais elle a peut-être entraîné 
une certaine paresse à aller au fond du problè- 
me. Nous esquisserons, en nous limitant toute- 
fois au domaine des faibles densités et à un type 
d'interaction, des généralisations de la formule 
de Balian et De Dominicis : le phénomène qu’ils 
ont découvert s'illustre davantage si on pro- 
cède à des rapprochements avec la théorie des 
phénomènes irréversibles. Les paramètres de 
l’équilibre thermodynamique disparaissent totale- 
ment de certaines formules : il est naturel que ces 
formules s’appliquent aussi bien hors de l’équi- 
libre qui ignore ces grandeurs qu’à l’équilibre 
même. C’est ce que nous espérons pouvoir mon- 
trer. 

Un exemple classique de rapprochements de 
ce genre est fourni par l’approximation de Kirk- 
wood. Imaginons sa généralisation : au lieu de 
l’écrire dans l’espace ordinaire, écrivons-la dans 
l’espace des phases. Elle est, intuitivement, aussi 
satisfaisante hors d’équilibre qu’à l’équilibre. On 
l'utilise effectivement dans la théorie des plasmas, 
qu’ils soient ou non à l’équilibre. 


IT. Relations algébriques. — Le problème exa- 
miné ici nécessite l’emploi de relations algébriques 
dont on verra qu’elles sont extrêmement simples. 
Elles développent un thème imaginé par Ursell. 
Nous suivons tout d’abord la présentation de 
Kahn et Uhlenbeck [6]. 

Faisant abstraction de tout problème de phy- 
sique, nous considérons un tableau de nombres. 
es N nombres à un indice, celui-ci allant de 

CAE 


UARUPAE .. .. AN.. 


For 
Puis 5 N(N — 1) nombres à deux indices pris 
dans la même suite allant de 4 à N : 


A2 A3 :.. Az ... Gij +. AN—1YN. 


AR 
Ensuite & N(N — 1) (N — 2) nombres à trois 
indices : 
Ci23 dira + ++ ik +. UN—2)(N—1)N- 
Et des suites avec un nombre d'indices de plus 
en plus élevé, jusqu’à la dernière suite qui ne 
comprend qu’un seul terme à A indices : 


d28...i...N: 


Dans chaque terme de ces suites, l’ordre des 
indices n’a pas d'importance. Ceci parce que nous 
posons par définition : 


or — 19; 
djik — Gijk, 
par exemple. Nous considérons 4,2 0t dy, die OÙ Gyix 


comme indistincts. 
Nous introduisons maintenant les combinaisons 


Ci, (1) 
Ciÿ — dij — di 4j, 
Cijk — Gijk — di Ax — dix Aj — Gjx Qi + 24 Gj a 
Ci = Gigi — Gije QU — Gij A — Gil j — ji di 
 dij O4 — ik A — y Ajk 
+ 2aij ax Q + 2x dj M + 2aiy a; a 
+ 24m di 4j + 2aj di dx + 24ajr di 4 — 64 4j ax A, 


dont on imagine aisément la généralisation pour 
un nombre quelconque d'indices. Les termes qui 
re comportent qu'un seul facteur algébrique ont 
le signe +, ceux qui en ont deux le signe —, 
ceux qui en ont trois le signe +, et ainsi de suite. 
Les facteurs numériques sont successivement : 


0 ! Au DA SRE 
On vérifiera alors aisément les formules sui- 
vantes : 
dis — Cid + Cire (2) 


d23 — C1 + C2 03 + C3 do + Ces 
G323a = C1 Gosa + Cyo Ga + Cia dos + Cia os 
+ Cios da + Cioa G3 + Cisa A2 + Ciosa 
Ggo...N = C1 G2s...N + 2 Ci 34... 
LE 22 Cios a5...N + -.. + Ciss...N. (3) 


Dans la dernière équation, la somme simple 
comporte N — 1 termes où le rôle de l’indice 2 est 
joué successivement par les indices 3,4, ... N. La 


1 
somme suivante comporte 5 N(N — 1) termes 
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où le rôle des indices 2 et 3 est joué successivement 
par tous les autres couples distincts que l’on peut 
prendre dans la suite 2, 3,... N. Et ainsi de suite. 

Ces formules, établies pour des nombres, restent 
valables s’il s’agit d'opérateurs, pourvu que 
ceux-ci commutent lorsqu'ils ne comprennent 
aucun indice en commun. Par exemple, nous sup- 
posons que 


dos Commute avec &s% 
mais nous n'avons pas besoin de supposer que 
12 commute avec Gug. 


Les définitions (1) et la relation (3) sont exploi- 
tées dans le paragraph? suivant. 


III. Développement de l’opérateur densité. — 
Le système étudié contient N particules identiques. 
Son hamiltonien s'écrit : 


HE D NS He De He OM 2 
RCI En me LU Do nel ie A 


Il comprend donc des termes relatifs à une parti- 
cule isolée et des termes d'interaction entre les 
particules prises deux à deux. La généralité de ce 
qui suit ne serait pas changée s’il existait aussi des 
termes liant les particules trois à trois, du type 


Wir 


ou quatre à quatre, par exemple, du type 
Wijr. 
Nous désignerons par /,,l’hamiltonien qu’aurait 


le système s’il ne comprenait que les particules 1 
et 2 


H55 —= H, + He + W 2. 
Nous définirons de même À ,,3 
Hsos = Hit HW Wis + Wos, etc. 


Le système est enfermé dans un thermostat ; son 
état est décrit par l’opérateur de Boltzmann : 


1 
Dx = Ni Ty B:2...N e—PX, (4) 


l, est un facteur de normalisation. B:2...xest 
la somme de toutes les permutations possibles 
des N particules, avec le signe + pour les permu- 
tations paires et le signe -L pour les permutations 
impaires, le signe + convenant au cas des bosons 
et le signe — au cas des fermions. On posera : 


CE 


Il est commode de pouvoir calculer B,,...x par 
récurrence : 
B:2...N (5) 
= (1 + ePys + Past... + eP5 +... + Pin) Bis... 
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P; étant l’opérateur de transposition qui permute 
seulement les particules à et 7. En outre, on note 
que : 


Di; 
Bi = 1 + ePis. 


Il est essentiel d'introduire explicitement B,2...n 
dans la définition de l’opérateur densité afin que 
la suite des fonctions orthonormées de référence 
ne comporte aucune condition de symétrie : toutes 
les conditions de symétrie sont contenues dans 
l’opérateur B. r. 

Posons maintenant : 


— 6H, LH 
A; —=e MM Ajj = e k Le 
—6H;;: —$H. 
A CM Ed ie à (6) 
di = Ai, Gij = Bi; Ai, 


aie = Bijx Aijx, ..., Qye...N = Bi2...N Au2...N (1) 

Quant aux grandeurs C, elles conservent leur 
définition. La formule (3) s'applique. Elle fournit 
le développement cherché de Dy. On introduit ce 
développement dans l’expression de l’opérateur 
densité simple : 


Uni = NTs...n DN. 
Il vient donc, compte tenu maintenant de 


l’équivalence des différents termes de chaque 
somme : 


NE 
Lx UN1 = C1 Tos...N Gas... N 


UN 1) TS Cris Ts ONG. NN 


1 
HS, AN — 2) Tes Cyos T'as... N Gas...N + ... 1 (8) 


Ce calcul est possible grâce à la nature des 
fonctions propres qui permet le calcul aisé des 
traces partielles. Par exemple : 


Tya...nN Ce Gsa...N = Cure Ta... N Asa...N- 


Le résultat (8), encore trop compliqué, se sim- 
plifie en supposant que maintenant le système est 
ouvert. Le thermostat est constitué par un immense 
réservoir de particules du type envisagé, qui 
communique avec le système étudié. La proba- 
bilité que ce dernier contienne W particules est 
maintenant : 


Ro (9) 


où z est l’activité et Æ un facteur de normalisation. 
Dans le système ouvert, l’opérateur densité à une 
particule est une moyenne qui s'obtient en pondé- 
rant par zx les opérateurs densités pm, : 


Us = ZIx Ua (10) 


No:7 
Ce calcul donne finalement l'expression cherchée : 
12 
by = 201 + 2 To Cie + Cie Ts Cios 


+34 Tosa Cuasa + +. (11) 


On obtient un développement en série entière 
par rapport à z. C’est la généralisation évidente 
du développement classique sur le même sujet. 


IV. L’opérateur densité double. — Il s’agit de 
l'opérateur u,, défini au mémoire [1]. Il nécessite 
un développement nouveau de l’opérateur D. La 
formule (6) faisait jouer un rôle particulier à la 
particule 4. Il faut maintenant mettre en évidence 
le couple 1-2. On pose : 


# 

C12 —= Gyo 
LA 

Ci23 — Aya — yo 3 $ (12) 
a 

Ci234 = A1234 — Uios da — Quog A3 — Uy2 A3a + 2012 A3 Ua 


De 


la même règle que précédemment s’appliquant sur 
les signes et les facteurs numériques et étant 
entendu que le couple 1-2 n’est jamais séparé entre 
deux facteurs distincts. Maintenant : 


# Je 
Be ON UC dei N DICO nt 


+ ZE Chase a56...n + -. (13) 

On conserve aux a les définitions (7). Le dévelop- 
pement (13) conduit, après les mêmes transfor- 
mations qu’au paragraphe précédent, à l'expression 
suivante de l’opérateur densité à deux particules : 


à LA ae 1 4 
Wie = 2? C2 + 2 T3 Ci23 + Sa TeinCissae ere (14) 


V. La suite du programme. — Nos raison- 
nements ont supposé successivement que le système 
étudié était un système thermodynamique fermé 
et plongé dans un thermostat puis un système 
thermodynamique ouvert, c’est-à-dire qu’alors le 
thermostat est capable de délivrer au système 
étudié, à travers une paroi semi-perméable, des 
particules de l’espèce désirée. Ces dispositions, 
passablement artificielles, sont commodes pour 
l'étude de l’équilibre thermodynamique. Ce que 
nous allons imaginer maintenant est tout aussi 
artificiel, mais mieux adapté à notre objet. Nous 
supposerons donc que le système est indéfini, le 
système s'étend à perte de vue. Même si les parti- 
cules s’agglomèrent quelque part, dans un léger 
creux de potentiel, pour former une phase dense, 
un liquide par exemple, on en trouve à toute dis- 
tance, plus ou moins raréfiées suivant les circons- 
tances. On ne précise pas si le nombre total des 
particules est fini ou infini : il pourrait arriver 
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encore que le milieu soit indéfiniment homogène. 
Le système constitue son propre thermostat. 

On conçoit que les formules (11) et (14) s’appli- 
quent encore au cas d’un tel système, mais les 
Opérateurs À,, A,o, A393, ... ont maintenant un 
spectre continu — non pas nécessairement tota- 
lement continu. L'intérêt d’une telle situation est 
que nous pourrons profiter maintenant des pro- 
priétés de l’opérateur de diffusion. 

Considérons d’abord un système ne contenant 
que deux particules. L'opérateur de diffusion est la 
limite suivante, si elle existe : 


2 ) exp. (ie ): (15) 


Cette limite n’existe pas pour n'importe quelle 
loi de force. Jauch et Zinnes [7] ont précisé les cas 
favorables à propos d’une simple force centrale. 
Nous nous plaçons ici dans le cas où une telle 
limite existe. Nous supposerons de plus que l’opéra- 
teur Q est unitaire : ceci a lieu si les opérateurs 
H, + H, et H,, ont le même spectre. Nous lais- 
serons entièrement de côté les généralisations que 
pourrait imposer le cas général. L'opérateur uni- 
taire Q,, transforme les fonctions propres de 
H, + H, en fonctions propres de 4, + H, +W,, 
en conservant la norme et l’énergie. [l en résulte 
la relation suivante : 


Q,, = lim exp. (- 1 
1=co 


exp [— B(Hi2)] = Que exp [— B(H: + H3 Q, (16) 
ou, avec nos notations : 


4e (04 A, 05. (17) 


En nous limitant aux termes en 2? — ce qui 
suppose qu’il s’agit d’un milieu peu dense, d’un 
gaz — il vient d’après (11), (14) et (17) : 


Wu = 241 + 2 To(By Q;: A: A0 (or — A; A:), (18) 


Hi2 — 3? Bis Quo A9 A1 Q%.. (19) 


A la même approximation, on résout la première 
équation par rapport à 74 : 


li T(Bis Oo Wii We On — Hi U), (20) 
et on élimine la même quantité dans la second» : 


Us — Bis Os Ua be On". (21) 

Dans la limite de l’approximation, les deux for- 
mules précédentes réalisent le programme proposé. 
L'activité et la température ne figurent qu’au 
premier membre de (20). Quant à la formule (21), 
qui exprime les corrélations, elle ne contient aucun 
paramètre thermodynamique. La généralisation de 
cette procédure ne comporte formellement aucune 
difficulté si on suppose qu’il existe des limites Q,:3, 
Q, 234, --. qui généralisent Q,,. 

Dans le calcul de la trace T°, 12 y a certainement 
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intérêt à couvrir le terme y, u, par l’opérateur B,, 
On écrira : 


TÂB19 Qi2 Ui On — Hi U2) 
= T3 Bio Qo Mi Ho Or! ae Mi U2) + T2 Piste (22) 


Mais 


Ta Pie ob = ui. 


On écrira donc 
2À1 = Wii — eui HET B\0 Hi Ha OR Se Ka Ua). (23) 


On comprend qu’en opérant de même avec les 
termes d'ordre supérieur, on puisse faire paraître 


Ua — oui + ui — eut + ..., 
c’est-à-dire le développement de 


Ma a Ua 
1 + ep 


C'est tout ce qu’on trouverait en l’absence d2 
toute interaction. On retrouve aussi correctement 
les formules relatives aux particules libres. Ce 
genre de transformations, qui assure probablement 
une convergence meilleure, fera l’objet d’études 
ultérieures. 


VI. Équilibre et régime variable. — Il est remar- 
quable que la formule (21) s'applique aussi en 
régime variable. Nous reproduisons ici pour la 
commodité l’essentiel d’une note sur ce sujet [8] : 
en régime variable, les opérateurs densités pu 
évoluent suivant les équations de Born et Green : 


x À 
ls di = Hits — bi Hi + TafWio bio — Wie Wao), (24) 


= His Vie — Wio Ho 


+ Tal(Wis + Wos) Los — Uioa(Wis + Wos)]. 


On suppose que les interactions sont à courte 
portée. Il en résulte, parce que la densité est 
modérée, que les traces T, ou T,, sont relativement 
peu importantes. Négligeons-les donc en première 
approximation. Posons : 


Ur exp(— 4), Une — exp U 


Il vient 
et) & Uiu:(0) Ur, 
Wo(t) © Us U2(0) Un. 


Nous formulons maintenant l’hypothèse ici fon- 
damentale qu’il n’y a aucune corrélation entre les 
particules à l’instant { — 0 ou du moins qu'il 
n’y a pas d'autre corrélation que celle imposée par 
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la nature des particules, bosons ou fermions. Cette 
hypothèse s’écrit : 
Uo(0) — B39 41(0) W(0). (25) 
Il nous est alors possible d'exprimer p,,(t) en 
fonction de u,(t) et de (tv) : 


bot) & Us Unis B9 Ut) Bolt) Us Us Une). 


Nous supposerons encore ici que l’opérateur 
Us Un Ur 


tend vers une limite quand le temps £ est suffi- 
samment grand. Le système doit être indéfini. 
Cette limite est Q. Notre résultat s'écrit : 


We = Bio Oo bo Qi. (26) 


Dans le cas d’un gaz peu dense, ce résultat est 
assez précis pour être incorporé dans le terme de 
rencontre de l’équation (1). On obtient une généra- 
lisation quantique de l’équation de Boltzmann : 


s Ù 
PAS Gi at 


Te Pal Oo io Of — Qi be QE Wie]. 
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La déduction repose sur une hypothèse de cha9s 
moléculaire : cette hypothèse est exprimée par 
l'équation (25). Ici cette hypothèse n’est demandée 
que pour l'instant initial, contrairement à ce qui 
a lieu dans le cas du traitement classique habituel 
de l’équation de Boltzmann. 

Les corrélations s’expriment donc par la même 
relation en régime d’équilibre thermique (for- 
mule (21)) et en régime variable (formule (26)) 
lorsque l’opérateur Q est unitaire puisque. alors : 


Sr (9 


Une question se pose alors : la formule statique 
(21) complétée par les termes d’ordre supérieur 
s’applique-t-elle encore aux régimes variables ? 
Cette généralisation, si elle était justifiée, pré- 
senterait un intérêt évident. Cette généralisation 
ne saurait s’appliquer, bien entendu, qu’à des sys- 
tèmes organisés depuis fort longtemps et à peu 
près abandonnés à eux-mêmes, c’est-à-dire peu per- 
turbés de l’extérieur. 

L'auteur remercie M. C. Bloch d’utiles remarques. 


Manuscrit reçu le 25 mars 1960. 
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L'EFFET PHOTOVOLTAIQUE DE SURFACE DANS LE SILICIUM 
ET SON APPLICATION À LA MESURE DE LA DURÉE DE VIE DES PORTEURS MINORITAIRES 


Par M1 AxDRÉE QUILLIET et M. Perer GOSAR (*) 
Laboratoire de Magnétisme et de Physique du Solide, à Bellevue. 


Résumé. — On a vérifié expérimentalement sur des échantillons de silicium une relation prévue 
par l’étude théorique : pour chaque longueur d’onde l'effet photovoltaïque de surface est propor- 
tionnel à ÆL/(1 + ÆL), où K représente le coefficient d’absorption correspondant, et L la longueur 
de diffusion des porteurs minoritaires dans l’échantillon. 

La relation étant vérifiée dans un petit intervalle de longueur d’onde dans les conditions expéri- 
mentales utilisées, on a pu déterminer des valeurs de la durée de vie des porteurs minoritaires en 
accord satisfaisant avec les résultats de la méthode de Valdes. 


Abstract. — An experimental verification of a relation predicted by theory was made on sili- 
con samples : ateach wavelength, the surface photovoltaic effect is proportional to ÆL/(1 + KL), 
where Æ is the corresponding absorption coefficient, and Z the diffusion length of minority carriers 


in the sample. 


. The relation being shown valid in a small wavelength range, values of the minority carrier 
lifetime could be determined, and these were in good agreement with the results given by the 


method of Valdes. 


Introduction. — Une différence de potentiel, +, 
existe généralement entre la surface d’un semi- 
conducteur et son intérieur (fig. 1), due à des 
niveaux d'énergie supplémentaires localisés à la 
surface. Lorsque l’on éclaire la surface du semi- 


SURFACE 
SEMICONDUCTEUR ' 


POTENTIEL 


Fig. 1. — Variation du potentiel 
au voisinage de la surface d’un semi-conducteur. 


conducteur, © diminue, et cette variation, appelée 
effet photovoltaïque de surface [1], [2], [3], dépend 
de la longueur d’onde de la lumière et s’exprime 
en fonction du coefficient d'absorption, de la durée 
de vie des porteurs minoritaires et de la vitesse de 
recombinaison en surface. 

Nous donnons ci-après une étude expérimentale 
du phénomène et son application à la détermination 
de la durée de vie dans le silicium. 


(*) Institut « J. Stefan », Ljubljana, Yougoslavie. 


PRINCIPE DE LA MESURE 


1° Expression de la variation de ©. — La den- 
sité de charge superficielle, o, est fonction de la 
différence de potentiel, +, du facteur de dopage, 
À — po/ni, (p, étant la concentration des porteurs 
minoritaires, #; celle des porteurs majoritaires, 
dans un échantillon de type », pour fixer les idées) 
et du facteur d’injection, « — Ap/nri, (Ap étant 
l’accroissement de la concentration des porteurs 
minoritaires dû à l’éclairement) : 


Gi ADN O0: 0) 


Lorsque la surface est éclairée pendant un temps 
assez court, on peut considérer la densité de charge 
superficielle comme constante, et la différence de 
potentiel de surface ®, prend une valeur @ telle 
que : 

F($: ÿ, O) — F(?, à, &). 


Si l’éclairage n’est pas trop intense, on peut 
admettre que + est une fonction linéaire de « : 


D—p) — Ko.x — k.Ap. (4) 


20 Mise en évidence de la variation de @. — La 
méthode de mesure utilise un contact capacitif avec 
l’échantillon (fig. 2). La différence de potentiel 
alternative apparaissant aux bornes du condensa- 
teur ainsi constitué et soumis à un éclairement 
périodique permet de déterminer @ — %. 

Pour que la lumière puisse atteindre la surface 
du semiconducteur, le diélectrique et l’électrode 
doivent être transparents. Cette dernière est une 
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plaque métallique perforée. Le schéma de montage 
est celui de la figure 2. L'ensemble est placé der- 
rière la fente de sortie d’un spectromètre Perkin- 
Elmer modèle 112 à simple faisceau et reçoit un 
signal lumineux modulé à 13 c/s. Une lampe élec- 


semiconducteur 
électrode 


F1G. 2. — Schéma du montage. 


tromètre E adapte l’impédance (À = 101% ohms 
à l’entrée, 20 000 ohms à la sortie). Le signal est 
reçu par l’amplificateur Perkin-Elmer modèle 81. 

L'emploi d’une plaque perforée a un inconvé- 
nient : des portions de l’échantillon ne sont pas 
éclairées et, en l’absence de diffusion, on observe- 
rait un effet très petit. Mais on peut obtenir un 
effet photovoltaïque important lorsque la distance 
et le diamètre des trous sont de l’ordre de la lon- 
gueur de diffusion. Cette méthode ne s'applique 
donc qu'aux cristaux de longue durée de vie. Il 
peut être utile d’augmenter la couche d’inversion 
par des traitements chimiques ; il faut éviter les 
surfaces polies mécaniquement et 1l est nécessaire 
de décaper chimiquement juste avant la mesure. 
Il est également utile de réduire au maximum 
l’épaisseur du diélectrique. 


30 Caleul de l’accroissement de concentration 
des porteurs minoritaires créés par la lumière. — 
La tension photovoltaïque de surface est, d’après 
l’équation (1), proportionnelle à l’augmentation 
totale de concentration des porteurs minoritaires 
au voisinage de la surface, Ap. La condition de 
continuité, exprimant la relation entre la diffu- 
sion, la recombinaison et la génération optique des 
porteurs minoritaires, est : 


2 

PP PE Ph)K.oe, (2) 
K est le coefficient d'absorption, p est l’augmen- 
tation de concentration des porteurs minoritaires 
à la distance x de la surface, D la constante de 
diffusion, + la durée de vie. P(à) est le nombre de 
paires électron-trou créées par la lumière de lon- 
gueur d'onde À. Le rendement quantique étant 
égal à un, P(A) est égal au nombre de photons 

absorbés. 
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T.S. Moss [4] a donné la solution de l’équation (2) 
en tenant compte de la recombinaison en surface, 
s, supposée faible. Dans le cas où l’épaisseur de 
l'échantillon est grande devant la longueur de 
diffusion, l'accroissement total de concentration 
des porteurs est : 


P()) TS KL a 
APT Der e le (3) 
(L: longueur de diffusion des porteurs minori- 
taires). 
Posons : 


LGV ARE 
VE Po CADET 


Ici À est une constante indépendante de À. 
On peut donc déterminer L par la formule : 


YIR=ETYEE A, (4) 
et par suite, T, puisque: 
x = L°]D. (5) 


4° Détermination du coefficient d’absorption. — 
Les valeurs de Æ sont déterminées sur une courbe 
d'absorption calculée point par point à partir des 
facteurs de transmission T et de réflexion À, mesu- 
rés sur un échantillon d'épaisseur connue d ; nous 
avons utilisé la formule approchée suivante : 


T = (1— R)°.exp (— Æd). 
PARTIE EXPÉRIMENTALE 


Des mesures ont été effectuées sur cinq pla- 
quettes de silicium monocristallin dont les carac- 
téristiques sont données dans le tableau 1. 


1° Détails opératoires. — La source était un 
élément Globar. 

Suivant l’amplitude (variable avec l’échantillon) 
du signal obtenu, on a utilisé diverses ouvertures 
de la fente d’entrée du spectromètre, correspondant 
à différentes largeurs de fente spectrale. Celles-ci 
sont données dans le tableau I. 

Deux plaques perforées en cuivre ont été utilisées. 
La plaque I, transparente à 36 %, portait des 
trous de diamètre moyen 0,15 cm, distants de 
0,125 em. La plaque IT, transparente à 25 ©}, 
portait des trous de diamètre moyen 0,045 em, 
distants de 0,065 em. Entre l’échantillon et la 
plaque perforée était insérée une feuille de mica 
d’une épaisseur de 10 microns. 


2° Variation spectrale de l’effet photovoltaïque 
de surface et du coefficient d’absorption. — Le 
signal photovoltaïque a été enregistré en fonction 
de la longueur d’onde. La figure 3 le représente 
sous forme de courbes normalisées pour un nombre 


NO, L'EFFET PHOTOVOLTAIQUE DE SURFACE 977 
TABLEAU I 

No Rér Tres ÉPAISSEUR RÉSISTIVITÉ DURÉE DE VIE DURÉ 

E Es s d(cm) p(Q@.cm) VALDES (us) ERSOnE A) SORA" 
1 a p 0,052 75 0,7 se 0,11 0,66 

2 FZ p 0,161 72,6 24,0 I 0,11 24,6 

3 DV4 p 0,035 0,5 3,5 Il 0,08 3,0 

[A DV6 n 0,034 0,4 _ Il 0,05 26,4 

5 » » » » » IT 0,03 27,8 

6 DV? n 0,052 0,4 7,8 Il 0,08 8,5 


de photons incidents constant. Des courbes d’ab- 
sorption ont été tracées pour les diverses ouver- 
tures de fente. Ces courbes diffèrent légèrement 


y [(unités arbitraires) 


—#— ——— 


—. + r _. — z 
0,90 092 0,94 0,96 0,98 100 1,02 1,04 ap} 


0 
082 084 085 088 


F1G. 3. — Courbes de réponse spectrale de l'effet photo- 
voltaïque de surface pour un nombrefde-photons inci- 
dents constant. 


entre elles et se rapprochent, pour les fentes les 
plus étroites, de la courbe publiée par Fan [5]. 


3° Détermination de la durée de vie. — A partir 
des données précédentes, on a construit point par 
point les courbes représentant y/XK en fonction de 
y (fig. 4). D’après l’équation (4), ces courbes 
doivent être des droites. De leur pente, on déduit 
L, puis la durée de vie + au moyen de la formule (6), 
en prenant D — 38 cm?/sec ou D = 13 cm?/sec 
pour le silicium de type p ou n, respectivement. 


4° Discussion des résultats et conclusion. — 
Nous constatons un bon accord entre la valeur 
donnée par la méthode de Valdes et celle que nous 
obtenons par l’analyse des courbes de la figure 4. 
La précision de nos mesures peut être estimée 
entre 10 % et 25 % selon la valeur de +. 

Toutefois, les portions linéaires des courbes de la 
figure 4 ne correspondent qu'à un domaine de 


longueur d’onde limité, dans lequel les signaux 


photovoltaïques mesurés sont en bon accord avec 


les prévisions théoriques. Pour les longueurs d’onde 


y (unités arbitraires) 
CARE] 2 3 L 5 6 7 
Fi. 4. — Pentes des droites y/X. 


plus grandes ou plus petites, la précision devient 
mauvaise et les résultats s’écartent notablement de 
la théorie : les valeurs expérimentales de y/X ne 
varient plus linéairement avec y (cf. par exemple 
courbe 3, figure 4). 

On améliorerait la précision de la mesure en 
choisissant le meilleur compromis possible entre 
deux exigences contradictoires : haute résolution 
et signal intense. 

Dans ce but il faudrait diminuer la fente 
d'entrée, utiliser une plaque à trous nombreux 
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et rapprochés, ou une couche métallique semi- 
transparente, et rendre éventuellement plus sen- 
sible le système d’amplification du signal. Des 


JOURNAL DE PHYSIQUE 


No 


durées de vie de l’ordre de 0,1 us sont certaine- 


ment mesurables. 
Manuscrit reçu le 11 avril 1960. 
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REVUE DES LIVRES 


ARZELIÈS (H.), avec la collaboration de MENDEz (R.), 
La dynamique relativiste et ses applications. II. Pro- 
blèmes de mouvement en dynamique du point faible- 
ment accéléré. Un vol. de 451 p., 16 X 25, avec 
187 fig. — Gauthiers-Villars, Paris, 1958, prix 60 NF. 
Les fascicules successifs des Etudes relativistes de H. Ar- 

zeliès sont appelés à faire autorité auprès de nombreux 
techniciens utilisateurs de la Théorie, à avoir aussi la 
meilleure influence sur les nombreux esprits réfractaires 
aux exposés abstraits mais désireux pourtant de com- 
prendre sans contre-sens. 

Le but général et le résumé du contenu de ce très beau 
livre sont excellemment présentés par l’auteur lui-même 
en les termes que voici : « Le présent volume désire fournir 
l’équivalent relativiste de nombreux ouvrages traitant de 
la dynamique du point. Il contient l’étude détaillée d’un 
certain nombre de mouvements qui se retrouvent fré- 
quemment dans les applications. Je réserve les accéléra- 
teurs de particules et l’optique corpusculaire pour des 
fascicules ultérieurs. De nombreux problèmes de mouve- 
ments relativistes ont fait l’objet d’articles, mais un exposé 
d’ensemble substantiel n’existe pas encore. Les utilisa- 
teurs ont cependant besoin d’un tel exposé. L'ouvrage 
s'adresse surtout aux physiciens et aux ingénieurs. Mais 
comme je me suis placé à un point de vue général, je traite 
parfois des questions dont je ne vois pas actuellement 
d'applications. Elles ont l’avantage de se ramener aux 
quadratures et constituent d’excellents exercices. » 

D'accord avec nous-même, H. Arzeliès considère que 
la définition générale de la force en Relativité est iso- 
morphe à celle qu’enseigne l’Électromagnétisme (force de 
Coulomb + force de Laplace). Ceci étant, il développe 
systématiquement des méthodes de calcul générales per- 
mettant de traiter les diverses classes de problèmes prati- 

> — 

quement utiles : champs E et B constants et uniformes, 

puis variables dans l’espace seul (avec nombreuses appli- 

cations : forces centrales et focalisation, rayons cosmiques 
et aurores boréales, effet corona, précession de Larmor. 


= 
oscillateur harmonique, magnétron); champ E dépendant 
du temps ou de la vitesse ; engins à réaction, astronau- 
tique ; chocs de particules relativistes (élastiques ou iné- 
lastiques) ; notion d’antiparticule selon Stueckelberg- 
Feynman. 

La proposition de M. H. Arzeliès relativement au forma- 
lisme 4-dimensionnel est connue, et parfaitement logique 
étant donné le but qu’il vise. Voici comment il la définit 
lui-même : «(le formalisme 4-dimensionnel) est très peu 
utilisé et j’en donne seulement deux exemples. Aucun 
technicien ne l’emploie. Très intéressant pour les questions 
épidémologiques et pour certaines théories générales, il 


n’est qu’une gêne en cinématique et en dynamique (appli- 
quées) ». 

Si, pour finir, nous louons la clarté du texte et des cal- 
calculs, l'abondance de diagrammes très soigneusemen7 
exécutés, le piquant d’une langue intéressante et vivante 
et le charme de préfaces « à la Bouasse » (mais de plus de 
portée philosophique), c’est afin d’inciter tout technicien 
electronique ou nucléaire, tout «honnête homme», à 
prendre connaissance de ce magistral ouvrage, en tout 
point digne des deux précédents. 


O. CoSTA DE BEAUREGARD. 


ARZELIÈS (H.), avec la collaboration de Henry (J.), 
Milieux conducteurs ou polarisables en mouvement. 
Un vol. de xLiv + 347 p. (16 X 25) avec 86 fig. — 
Gauthier-Villars, Paris, 1959, prix 58 NF. 

Ce quatrième volume de l’encyclopédie relativiste de 
H. Arzeliès est le très digne successeur des précédents ; il 
se trouve être aussi le successeur de l’électrodynamique des 
corps en mouvement de H. Minkowski. 

M. H. Arzeliès montre à plusieurs reprises que les pro- 
blèmes d’électrodynamique caractérisés par son titre ne 
peuvent même pas être correctement posés sans un recours 
explicite aux principes de la Relativité restreinte (voir par 
ex. p. 88-89) ou de la Relativité générale (voir la 2€ partie). 
Il montre aussi que, contrairement au sentiment répandu, 
la théorie correcte, et toujours subtile, des effets Cullwick, 
Rowland, Crémieu-Pender, Wilson, Rœntgen-Eichenwald, 
n’était pas toujours acquise avant sa propre intervention. 
Il attire l’attention sur le fait que les nouveaux matériaux 
magnétiquement ou électriquement polarisables dont on 
dispose aujourd’hui permettraient de pousser ces expérien- 
ces macroscopiques beaucoup plus loin qu’on ne l’a fait 
encore. En bref, et contrairement à ce qu’on croirait à 
première vue, il reste là des recherches à faire soit théo- 
riques, soit expérimentales. 

Et l’on découvre, en suivant les raisonnements et les 
calculs de M. H. Arzeliès, que l’électrodynamique des corps 
en mouvement constitue une mine quasi inépuisable de 
problèmes de deuxième ou de troisième cycles. Parmi 
ceux qu’il s’est attaché à traiter, citons au hasard celui de 
aimant sphérique ou cylindrique en rotation uniforme. 

Plusieurs appendices mathématiques (calcul tensoriel, 
fonctions spéciales) ou physiques (grandeurs et unités) 
ainsi que des index d’auteurs et de notations complètent 
cet ouvrage magistral. 

Comme dans ses trois précédents ouvrages, l’auteur 
donne une bibliographie abondante des sujets traités. 
L’on reste émerveillé de l'étendue de son ‘information de 
premiere main, tout autant que de la lucide rigueur avec 
laquelle il fouille les problèmes jusqu’au tuf. 


O. CosTA DE BEAUREGARD. 
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NORMALISATION UNIFORME DES ONDES PLANES DE SPIN 1 /2 
SOUS LES OPÉRATIONS D'’INVERSION 


Rare ET/OU: 


Institut Henri-Poincaré, Paris. 


Résumé. — On examine la possibilité d’une normalisation uniforme des ondes planes mono- 
chromatiques par rapport aux opérations d’inversion. On démontre que de telle normalisation est 
réalisable pour chacune des huit inversions unitaires, mais non pour les inversions antiunitaires 
sauf Cet PT M. Une normalisation uniforme est même possible pour plusieurs inversions simulta- 
nément, en particulier pour les quatre opérations CPT, M, CT, PM. Grâce à l’existence de norma- 
lisation uniforme des ondes planes pour les inversions unitaires, on trouve enfin des opérateurs du 
champ spinoriel qui sont invariants par rapport à ces inversions, notamment à CPT, M, CT, PM. 


Abstract. — The possibility of a uniform normalization of monochromatic plane waves is exa- 
mined with respect to inversion operations. Itis shown that such a normalization is realizable for 
each of the eight unitary inversions, but not for the anti-unitary ones, except C and PTM. A uni- 
form normalization is even possible for many inversions simultaneously, in particular for the four 
operations CPT, M, CT, PM. Owing to the existence of uniform normalization of plane waves 
with respect to unitary inversions, we find finally spinor field operators which are invariant under 


TOME 21, JUILLET 1960, paAce 579. 


these inversions, especially under CPT, M, CT, PM. 


1. Introduction. — Le champ de spin 1/2 étant 
une théorie fondamentale en physique des parti- 
cules élémentaires, nous nous sommes beaucoup 
intéressés à sa structure interne vis-à-vis du groupe 
de Lorentz Let plus particultèrement du groupe de 
symétrie composé d’inversions de la charge C, de 
l’espace P, du temps T'et dela masse M, groupe d’in- 
version à seize éléments (opérations d’inversion) : 


MP PT, CG, CP, OT, CPT, 
M, PM, TM, PTM, CM, CPM, CTM, CPTM }. 


Dans des articles précédents [1], [2], [3], nous 
avons établi, entre états propres d’une particule 
de spin 1/2, toutes les relations engendrées par les 
opérations de ces groupes, en particulier par celles 
du groupe d’inversion. À ce sujet, nous avons étudié 
comme exemples les ondes planes monochroma- 
tiques, les ondes sphériques et le problème de force 
centrale. 

Dans ces études antérieures, nous n’avons consi- 
déré, entre états propres, que les relations quali- 
tatives sous les opérations L, C, P, T, M. Ces rela- 
tions n’ont été précisées qu’en ce qui concerne les 
constantes de mouvement, soit les nombres quan- 
tiques des états propres considérés. Nous n’avons 
pas donné des relations précises des fonctions 
propres elles-mêmes, qui n’ont été identifiées sous 
les inversions qu'aux constantes de normalisation 
près. 
En relation avec une étude sur la forme parti- 
culière V — A des interactions faibles des fer- 
mions [4], nous avons proposé des interactions 
apparemment non locales qui se réalisent par 
«images d’inversion » par rapport à certaines opé- 


rations d’inversion, notamment CPT et M [5] ou 
PM et CT [6]. À ce sujet, nous devons bien préci- 
ser les relations exactes entre fonctions propres 
sous l’opération de ces inversions. Nous nous inté- 
ressons surtout à la possibilité d’une telle relation 
exacte qui, pour une inversion donnée, est indépen- 
dante des états propres considérés, donc indépen- 
dante des nombres quantiques. Si les fonctions pro- 
pres sont normées de cette manière par rapport à une 
certaine opération, nous appellerons une « norma- 
lisation uniforme » de ces fonctions propres vis-à- 
vis de cette opération. Évidemment, si les fonc- 
tions propres sont normées de façon «uniforme » 
pour une certaine inversion, il ne sera en général 
plus possible qu’elles le soient encore pour d’autres 
inversions. 

Dans ce qui suit, nous examinerons la question 
de normalisation uniforme dans le cas des ondes 
planes monochromatiques. Nous démontrerons 
l’existence de telle normalisation en cas d’inver- 
sion CPT, M, CT ou PM, résultat que nous avons 
supposé et dont nous nous sommes servis déjà dans 
des études antérieures [5], [6]. 


2. Définition des ondes planes et classification 
des inversions. — Nous définissons les ondes planes 
monochromatiques par (avec À — c — 1) 


PH) (6) (p ur t, m) — u)(s){p, m) ei(p.rTEt) (1a) 
PH) (p ; »,t, — m) = up, — mjeilprFE),  (1b) 
où E = +V/p? +m?,et s — + 1, valeurs propres 


de l’hélicité o, — ©6.p/p. 
Puisque nous allons considérer des relations 
exactes entre fonctions propres, il faut bien préci- 
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ser les constantes de normalisation des fonctions u. 
D'abord, nous les supposons normées de la manière 
usuelle, 

— 6, &) 00,10). (2a) 
utE)(s)f(p, — m) uÆ)(6)(p’,— m) — Ô(s, s’) Ô(p,p'). (26) 


u(Æ)(s)f(p, m) uGE)(s)(p", m) 


Même ainsi normées, les w ne sont encore déter- 
minées qu’à un facteur de phase près. Nous repré- 
sentons par uf (p, + m) les fonctions propres nor- 
mées selon (2) dont les facteurs de phase sont bien 
définis d’une manière quelconque, de sorte seule- 
ment que soient vérifiées les identités (1) 


US) (p, + m) — RURIERE P, F mn). (3) 
Dans ce cas, nous pouvons écrire 
(+ 
uÿ (P, + m). (4) 
Les € sont des facteurs de phase arbitraires 


CH) *(p, + m) CHG)(p, + m) = 1. (5) 


uH(p, + m) = CÆS(p, + m) 


En prenant une représentation particulière des 
matrices de Dirac, nous donnerons les expressions 
explicites des fonctions w, dans l’Appendice à la 
fin de l’article. 

Selon leur façon d’application, les seize éléments 
du groupe d’inversion se classent d’abord en huit 
opérations unitaires (7, P, CT, CPT, M, PM, 
CTM, CPTM) et huit opérations antiunitaires 
CHPPTSEC CP ETIENNE PECUÉeCPM) Pour 
faciliter l’étude ultérieure, nous divisons chacune 
des deux classes en opérations indépendantes de 
M et opérations comportant M. Nous aurons donc 
quatre ensembles d’inversions : 


= UR LNCARNONAIE (6a) 
AM = M, PM, CTM, CPTM, (6b) 
B=HRNPIACANCE (6c) 
BM = TM, PTM, CM, CPM. (6d) 


Conformément aux définitions, faites dans des 
articles antérieurs [1], [2], des états propres engen- 
drés par les inversions, nous pouvons définir ici, 
en cas d’ondes planes, 


(À) u(m) = l'4 u(m), (7a) 

(AM) u(m) = T'4m u(— m), (7b) 

(B) u(m) = T8 uT(m) — Qp u*(m), (7c) 

_ (BM)ulm) = Tam uT(— m) = Qpmu*(— m), (7) 
ou 

O3 = ls v?, Ogu = leu y. (8) 


L'indice supérieur Ÿ désigne la transposition des 
matrices. Les lpeuvent être des matrices unitaires 


(1?) En effet, u(Æ)(s (p,+m) et uP(9 (— p, + m) 
obéissent à des équations identiques, donc 
u(Æ) (s) (D, 2 m) & UF) (5) (— P,F m), 
# désignant l’identité à constante de normalisation près. 
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et her.aitiennes. Dans ce qui suit, nous les suppo- 
serons toujours ainsi : 


DT (9) 
Les Q sont alors unitaires : 


Q—1 = Qt. (10) 


Dans les paragraphes suivants, nous examinerons 
si des normalisations uniformes des ondes planes (1) 
sont possibles pour certaines opérations d’mver- 
sion et, dans l’affirmative, nous montrerons com- 
ment se réalisent ces normalisations. 


3. Normalisation uniforme par rapport aux 
inversions A.-— D’après (7a), (3) et (4), nous avons | 


(A) u(m) = T4 um) = aa(m) Cm) CAY- Um) u(4)(m). (11) 


Ici u'4) (m) représente l’état propre engendré de 
u(m) par l’inversion À [1]. Si les l, prennent 
comme représentation de simples matrices usuelles 
à éléments + 1 ou +1, les coefficients a4 (m) 
seront alors aussi de ces nombres simples, + 1 ou 
+ à, En tout cas, a (m) as (m) = 1. De l’équa- 
tion (11), on tire, à l’aide aussi de (9), 


(À) utA4){m) = T'4 u(A)(m) — aa (m) Am) C—1(m) u(m). 
(12) 
Nous pouvons alors mettre, sans incompatibi- 
lité quelconque, 
aa(m) Une) LAYU(m) = 44 (m) CA (m) Cm) = am), 
où X7 (m) Ka (m) = 1, X4 (m) étant indépendant 
des nombres quantiques de w(m) et uw (m). Il 
est d’ailleurs évident que 


xa(m) = + 1 ou —1. 


Nous aurons dans ce cas 
(A) u(m) = (+ ou —) ut4)(m), 
(A4) ut4(m) = (+ ou —) u(m). 


(13) 


D'une manière semblable, nous pouvons avoir 
aussi 
(A}u nm) = (ou) ut se) 


Au) (13) 


En conclusion, les ondes planes monochroma- 
tiques peuvent être normées de façon uniforme 
dans le cas des inversions À, et la normalisation 
uniforme se réalise, en cas de chaque opération, par 
seulement deux possibilités comme celles indiquées 
dans (13) et (14). Les facteurs de phase € doivent 
vérifier correspondamment les relations 


a am) fm) A Um) = + 1 où —1, (15) 
aa(— m) C(— m) CA m) = + 1 où — 1. (16) 


En prenant la représentation particulière (A-1) 
pour les matrices de Dirac (voir Appendice), nous 
donnons ci-après les normalisations uniformes des 
ondes planes sous les opérations d’inversion 1, P3 
CTACPT 


N°7 
TI. La normalisation est toujours uniforme. 
P. Nous avons l'? — y4. Alors 
(P) uÆs)(p, m) — 
+ CH6)(p, m) CHO—1(— p, m) WHO (— p, m), 
: (P)uÆG(p, — m) — 
+ CHE (p, — m) CHI p, — m) up, — mm). 


La normalisation uniforme sera possible, si les 
€ vérifient les relations 
e CH (6)(p, m) CH 6)—-1(— p, HO = SRE (17a) 


2e CÆ)(s)(p, — m) CŒH(S—1(— p, — m) = + 1 ou —1. 
(176) 


On aura alors correspondamment 
(P) up, m) = (+, ou —) uw p, m), 
(P) ut2)(p, 


(18a) 
} uÆ)(9)(— p, — m). (18b) 


m) = (+ ou 
MCE = 1Y): Y.. Ainsi 
(CT) u#6)(p, m) — 
+ (DO (p, m) CPE—U(p, m) uPW{(p, m), 
(CT) uE)6)(p, — m) — 
+ 1 CDO(p, — m) CP—1(p, — m) UPS (p, — m). 


La normalisation uniforme sera donc possible, 
sl 


+ it (p, m) CPO—Up, m) = +1 ou —1,  (19a) 
> + c CE)6)(p, — m) CME—-1(p, — m) — + 1 ou —1. (19b) 
Alors correspondamment 

(CT) u)S)(p, m) = (+ ou —) uF6)(p, m), (20a) 
» (CT) u&)®(p, — m) = (+ où —) FO (p, — m). (206) 

CPT. Ona ici l'or = Ys- Alors 
(CPT) up, m) 

= EH (p, m) (POI p, m) uFO(— p, m) 

(CPT) ut) )(p, — m) 

= LH6)(p, — m) CPO—U(— p, — m) F6) (— p, — m). 


La normalisation uniforme sera réalisée, si seu- 
lement 


CÆ)6)(p, m) CMPG—-I(— p, m) = + 1 ou —1, (21 a) 

CÆNS)(p, —m) CPG—-1(— p, — m) = + 1 ou — 1. (21b) 
On aura dans ce cas 

(CPT) us) (p, m) = (+ ou —) uPW(— p, m), (224) 


(CPT) us (p, —m) = (+ ou—) uF\9(—p,—m). (22b) 


Toutes les formules citées ci-dessus restent 
valables quand p — —p. 


4. Normalisation uniforme par rapport aux 
inversions AM. — Le problème est pareil à celui 
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traité dans le cas des A. En fait, d’après (7b), 
(3) et (4), nous avons 
(AM) u(m) = T4m u(— m) 

= GAmÛm) — m) CAM) uAM (1m), (23) 
puis évidemment 


(AM) uAM)(— mn) = T'4m ut AM)(m) 
= aam(m) CAM) (1m) Cl— mm) u(— m). (24) 
Ici 4 (m) est l’état propre engendré de 
u(m) par AM [2], et a%m (m) aam (m) = 1. Tout 
comme dans le cas des À, on peut mettre ici, sans 
incompatibilité quelconque, 
aam(m) C(— m) CAM-Um) 
— daim) CAM (m) C4 m) = + 1 où — 1. 
Ainsi, dans le cas des À M, la normalisation uni- 
forme des ondes planes est possible. Elle se réalise 


d’ailleurs dans deux façons différentes : Quand les 
facteurs de phase € vérifient les conditions 


aam(m) C(— m) CAM)-I{m) = + 1 ou — 1, (25a) 
aam(— m) Cm) CAMI-I— mm) = + 1 ou —1, (25b) 
nous aurons les résultats correspondants 
(AM) u(m) = (+ où —) wt4Um)] pen 
(AM) uAM){m) = (+ ou —) u(m), 
Rd 
(AM) u4M)(— m) = (+ ou —) u(— m). E 


Avec les matrices (4-1), nous pouvons démon- 
trer que la normalisation uniforme des ondes planes 
se réalise par rapport aux inversions M, PM,CTM, 
CPT M de la façon suivante. 

M. Nous avons l'y = y;. Alors 


(MjuÆ)6)(p, M) 
— C(Æ)(s)(p, —— mn) CH) —1(p, m) u(Æ)(s)(p, m), 


(M)uH)6)(p, = m) 
—— C(Æ)(s)(p, m) CH))—-1(p, == m) u(Æ)6)(p, —— m). 
La normalisation uniforme est possible ici, si 
CH) 6)(p, — m) (HG Up, m) = + 1 ou — 1. (27) 
On a alors 
(Mju®(p, m) = (+ où —) # HO (p, m), 
(Mjut(s)(p, — m) = (+ ou —) uÆ) 5 (p, — m). .(28b) 
PM. Ici, L'py = VYa Ys. Ainsi 
(PM)uÆ)6)(p, m) — 
+ DO (p, — m) CSI p, m) HO (— p, m), 
(PMj}u&)9(p, — m) — 
me LCÆ)(S)(p, m) CHN(S)—1(— P, + Im) uÆ)(s)(— p, — m): 


(28a) 


Si nous mettons 
Li CÆ6(p, — m) CHI p, m) = +1 ou — 1, (29) 
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nous aurons la normalisation uniforme, 


(PMjuÆ)s)(p, m) — (+ ou —) uÆO(— p, m), (30a) 
PMjuÆ0(p, — m) = (4 où —) dE (— p, — m). 
(30b) 


CTM. On a 101 Lorm = et puis 
(CT M)uÆ(s)(p, m) — 


(CT M)uÆ)6)(p, = m) — 
+ CÆ)(s)(p, m) CS) 1(p, — mj)utP)(6)(p, — mn). 


Nous mettons 


+ CE) (p, — m) CP—U(p, m) = + 1 ou —1, (31) 
et nous aurons 
(CT MjuÆ)6)(p, m) = (+ ou —) uF)6)(p, m), (324) 
(CT MjuÆ)s)(p, — m) = (+ ou —) uF)6)(p, — m). (32b) 
CPTM. Dans ce cas, Poppy = 1. Donc 
(CPTMj)uÆ)s)(p, m) 
= CH 6)(p, — m) (PO—1(— p, m) F9 (—_p, m), 
(CPTM)u#®)(p, — m) 
= CH)(s)(p, m) CH (s)—— D, ITS m) u(F)(s) (— D, —— Im). 
La normalisation uniforme se réalise, si l’on 
met 
CH S)(p, — m) CPE p, m) = + 1 ou — 1. (33) 
Nous avons alors 
(CPTM)juÆ(S (p, m) = (+ ou —)u PS {(— p, m), (34a) 
(CPTM)uÆSs)(p, — m) = (+ ou —) PO (— p, — 7m). 
(34b) 


Tous les résultats restent valables aussi pour 
st Lo 


5. Normalisation uniforme par rapport aux 
inversions B. — Le problème 1c1 est un peu diffé- 
rent de celui des À ou À M, car nous avons à traiter 
des opérations antiunitaires. Selon (7c), (3) et (4), 
nous avons 
(Bju(m) = Qg u*(m) — bg(m) C*(m) CB Um) utB)(m), 

(35) 
où 48 (m) est l’état propre engendré de u(m) par 
BH] et 05m) 06, (m) = 1."De’(35)/ontre, à 
l’aide de (10), 


Q$ utB)*{m) — BR 


(m) CB)*(7n) tm) u(rmn). 


Ici deux cas différents se présentent suivant que 
la matrice Q3 est symétrique ou antisymétrique. 

19 O3 symétrique, Q3 — Qf. L’équation (36) 
devient 


Bu B){m) = Qz ut B\*{m) 


(36) 


= bg ‘(m) CB)#{m») Cm) u(m).  (37a) 


+ (HO (p, — m) CPS—-IU(p, mu(P6(p, m), 


On peut mettre évidemment 
bglm) C*{m) Em) = es, 


os étant une phase réelle et arbitraire indépen- 
dant des nombres quantiques. Alors (35) et (37a) 
s’écrivent 


(38a) 


(Bu(m) — eieB u(B)(m), 


(BjutB(m) — 588 u(m). CE 


On peut avoir également 
be men er) rer, 


avec le même wz que us (38a). Nous aurons 
alors aussi 


(38b) 


(Bju(— m) — eieB wB(— m), 


(BjutB)(— m) — ei9B u(— m). (39) 


Ainsi une normalisation uniforme des ondes 
planes est possible pour chacune des opérations B, 
si Q3 est symétrique,.la phase 3 étant d’ailleurs 
arbitraire. 

20 Q, antisymétrique, Qz = — Of. 
cas, au lieu de (37a), nous avons 


(ButB){m) = Qp u(B)*(m) 


Dans ce 


= — bg “(m) CB)* (mn) Clim) um).  (37b) 


La normalisation uniforme des ondes planes n’est 
pas réalisable ici. En mettant (384, b), nous avons 
en effet 


(Bjulm) — eie8 ut(B)(rn), 

(B) utB){m) — — eisB u(m), QE 
(Bju(— m) = eiv8 uw B)(— m), 
(ButB)(— mn) = — eivB u(— m). | F9 


En prenant les matrices (A-1), on voit que 
Q, est symétrique, tandis que Qy, OQpr et Qop 
sont antisymétriques. Nous résumons ci-dessous 
les normalisations des ondes planes dans chacun 
de ces cas. 


T. Nous. avons Qi), y» —— Or, anti- 
symétrique. La normalisation uniforme n’est pas 
réalisable. En effet, on a (avec $s — — s) 
(TjuÆ9(p, m) 

— Si CÆ)(s)*(p, m) CH) (s—1(— P, mm) u(Æ)(s)(— P, m), 
(T}uÆs)(p, — m) — si CES *#(p, — 7m) 
CGI p, — m) u EE (— p, — m). 


Ici, on ne peut poser que des conditions comme 


CES) *(p, m) CHO—1(— p, m) = eier, (41a) 
LCHUSE(p, — m) CHE p, — m) = eisT. (41b) 
On aura donc seulement 
(Tu (p, m) = s eier uH@(— p, m), (424) 
(Thu (p, — m) = s eier uEHE@(— p, — 7 (42b) 


La normalisation dépend ainsi de l’hélicité. 


NO 
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PT. On à Qpr = is Ya = — Or, antis ymé- 
trique. 


(PT)ut(s(p, m) = 


£ SUD (p, m) LHG—1p, m) dE E(p, m), 


1" (PT)u (Sp, — m) — 


F Si CEUD*(p, — m) CH O—1(p, — m) HO (p, — m). 


On peut seulement mettre 


+ 1CŒO*(p, m) CHG—UWp, m) = eivpr, (43a) 
LC (p, — m) HG Up, — m) — eiepr.  (48b) 
Par conséquent, 

(PT)uÆ(s)(p, m) — s eiepr u +6) (p, m), (44a) 
(PThu)s(p, — m) = seiepr uÆ\9(p, — m).  (4kb) 


La normalisation dépend également de l’hélicité. 

CDans ce, cas, Qj — 1iy, — 0, symétrique. 
La normalisation uniforme est donc réalisable. En 
effet, nous avons 


(C}u)(s)(p, m) —= 
+ si CH(S)*(p, m) PE I p, muPO(— p,m) 
(CjuÆ(s)(p, — m) = + si CES *(p, — 1m) 
CG (— p,— m) uPE(— p, — m). 


Nous pouvons poser ici, sans incompatibilité, Les 
conditions 


£ si HO%(p, m) CPE p, m) — citc, 
F si CEE (p, So m) CRETE PS m) De CVS. 


Nous aurons alors la normalisation uniforme, 


(Chutt(s(p, m) = eirc uPH(— p, m), (464) 
(CjuÆ(s)(p, — m) = eiec uFIG(— p, — m). (46b) 
CP. On a ici Qop — Ye Ya — — fr, antisymé- 


trique. 
(CP)uÆ)(s)(p, m) 
= s (H(9*(p, m) CPE Ip, m) uPG(p, m), 
(CP)uæ&)(s)(p, — jf) 
= s CES #(p, — m) (PG—U(p, — m) u PH (p, — m). 


Évidemment, on ne peut poser ici que des 


| conditions comme 


CHS*(p, m) CPE —U(p, m) = eivcr, (474) 
CS #(p, — m) CPE —U(p, — m) — eifcr.  (4k7b) 
Nous aurons donc seulement 
(CP)uÆ®)(p, m) = s eircp u(FS(p, m), (48a) 
(CP)u&(s)(p, — m) = s ei9cp u(FI@(p, — m). (48b) 


La normalisation dépend encore de l’hélicité 
comme dans les cas T et PT. 
Tous les résultats restent valables aussi pour 


Pp — — P. 
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. 6. Normalisation uniforme par rapport aux 
inversions BM. — Le problème ici est pareil au 


cas des inversions B. Selon (7d), (3) et (4), nous 
avons 


(BMju(m) = Qgm u*(— m) 
— bpm(m) CX(— m) EMI) u(BM)(m), (49) 


où 4% (m) est l’état propre engendré de w(m) 
par BM [2], et bé (m) by (m) — 1. De (49), on 
tre, à l’aide de (10) 


QT, u(BM*{(m) — DE (mn) CBMIE (mn) CI — mm) u(— m). 
(50) 


Comme dans le cas des B, deux cas différents se 
présentent. 


19 Qyx symétrique. Nous avons alors 
(BMjuBM(— mm) = Qgm uBM*(m) 


= Em) EME) CU — m) u(— m). 


? 


(514) 


La normalisation uniforme est réalisable dans 
ce cas, car on peut bien mettre 


bpmlm) É(— m) CBM Um) — eirpm, (52) 


avec une phase réelle ozy indépendante des nom- 
bres quantiques. Puis, on a 


(BMulm) = eirpm u BM)(m), 


(BM)uBM(— yn) = eiepmy u(— m). 


(534) 
(53b) 

29 Oz antisymétrique. Au lieu de (514), on a : 
(BMjuBM(— m) = Qpy uBM*(m) 


a 0 n) CM rm) Sn) ur) (51) 


Une normalisation uniforme est donc irréali- 
sable, car, avec (52), on a seulement 
(BMu(m) = ei®8m u(B8M)(m), 


(BM)u(BM(— m) — — eiPpm u(— m). 


(544) 
(54b) 


Mais cette normalisation n’est différente que 
lorsque m —>— m. Elle ne dépend pas des nombres 
quantiques. C’est une normalisation « quasi-uni- 
forme ». 

Avec (4-1) pour les matrices de Dirac, nous 
donnons ci-dessous les schémas de normalisation 
des ondes planes pour les cas TM, PTM, CM, 
CRM 

TM. Nous avons Qypy = Yo Ya = — Om, anti- 
symétrique. 


(TM)uÆ\s)(p, m) 
— 5 CH *(p, — m) CHI p, m) uMO(— p, m). 


(TM)uÆ)(s)(p, — m) =. & CŒÆ(S)#(p, m) 
CH) — p, = Im) u(Æ)(s) (— D, — m) 
Nous mettons 


s (HS *#(p, — m) CHE —1(— p, m) = eivrm, (55) 
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et nous avons 

(TMjuH9(p, m) — eierm uHO(— p, m), (564) 

(TMjhuD®(p, — m) = — eirm uÆG(— p, — m). (56b) 
C’est une normalisation quasi-uniforme. 
PTM. Ici Qpru = 2Ye = Qfry, symétrique. 

On a 

(PT Mu 6) (p, m) — 

ENST CH)(s)*(p, = Im) CÆS)—1(p, m) uCE)(s)(p, 70) » 

(PTM)uÆ(s) (p, — m) — 

si CEO#(p, m) CHE—U(p, — m) HE (p, — m) 


On peut mettre,sans incompatibilité quelconque, 


+ si CES %(p, — m) CHE Ip, m) = eivprm, (57) 


et la normalisation uniforme est réalisable de la 
manière suivante, 


(PT M)uÆ\(s'(p, m) = eitpry us) (p, m), (58a) 
(PT M)uE)s(p, — m) = eifprm uÆ(s)(p, — m). (58b) 
CM. Dans ce cas, Qom = iy5 Ye = — Qu, 
antisymétrique. 
(CM)uÆ)(s)(p, m) — 
Æ SL CS) *(p, — m) CPG—I— p, m) u(F\s) (— P, m), 


(CM)u+s)(p, == m) = 
+ Si CH *(p, m) CPI — p, = m) us (— p, = m). 


En mettant 


+ si CH) *(p, 2m) CPI p, m) = eivcm, (59) 
on a une normalisation quasi-uniforme, 
(CM)uÆ)(s)(p, m) — eiPom uPS)(— p, m), (60a) 
(CM)u(S)(p, — m) = — eirom PS (— p, — m). (60b) 

CPM. Nous avonsici Qopy = lY1 Ys = — Qpy, 


antisymétrique. 
(CPM)uÆ)5(p, m) 
= si (HO *(p, — m) (PO—1(p, m) uPE(p, m), 
(CPMju#9(p, — m) 
= si CHO*(p, m) CPE—U(p, — m) u PE (p, — m). 
Si nous mettons 
si CES *#(p, — m) (PO Ip, m) = eircpm, (61) 
nous aurons une normalisation quasi-uniforme, 
(CPM)uÆs)(p, m) = eircpm uPS)(p, m), (62a) 
(CPMjuÆs)(p, — m) = — eivopm u(F){p, — m). (62b) 
ci-dessus restent 


Tous les résultats donnés 


valables pour p — — p. 


7. Normalisation uniforme simultanée pour plu- 
sieurs opérations d’inversion. — En général, si 
une normalisation uniforme est réalisée par rapport 
à une certaine inversion, elle n’est plus réalisable 


No 


pour d’autres inversions. Il est pourtant possible 
qu’une normalisation uniforme soit réalisable pour 
certaines opérations d’inversion simultanément. 
Nous examinons maintenant cette question. 

Nous nous bornons ici aux opérations unitaires. 
Par une inspection des conditions (19) et (21), une 
normalisation uniforme est possibl: simultanément 
pour CT et CPT. Elle n’est pas possible pour P 
et TR oU CRT 

On démontre également qu’une normalisation 
uniforme est réalisable simultanément pour M et 
PM et aussi pour CT M et CPTM, mais non pour 
M et CTM, etc. Ç 

En considération des interactions à «images 
d’inversion » [5], [6], il nous intéresse surtout 
d’examiner si une normalisation uniforme est simul- 
tanément réalisable pour CPT et M ou pour PM 
et CT. En fait, selon (21) et (27), une normalisation 
uniforme simultanée pour CPT et M entraine les 
conditions ) 


CH(9(p, m) = (+ où —) PO (— p, m) 


= (+ ou —) CE(s(p, — m), (63a) 
CH (p, — m) = (+ ou —) {P9(— p, — m) 
= (+ ou —) Es (p, m). (636) 


Ces conditions sont bien compatibles, et des 
facteurs de phase € existent donc pour qu’une nor- 
malisation uniforme soit réalisée pour CPT selon : 
(22) et à la fois pour M selon (28). 

D'autre part, d’après (19) et (29), on démontre 
également qu’une normalisation uniforme simul- 
tanée est réalisable pour CT selon (20) et pour 
PM selon (30), car les conditions 


LH (p, m) = (+ où —) à iCPE(p, m) 
— (+ ou —) + 1 CÆS{— p, — mm), (64a) 
Ep n) ou re pr) 
= (+ ou —) + : CÆ(S(— p, m), (64b) 


sont bien des relations compatibles aussi. 

Il existe même des facteurs de phase € qui vér.- 
fient à la fois les conditions (63) et les condi- 
tions (64). Une normalisation uniforme est done 
réalisable simultanément pour les quatre inver- 
sions CPT, M, CT, PM selon les schémas (22), 
(28), (20), (30) respectivement. 


8. Normalisation des opérateurs de champ par 
rapport aux inversions unitaires. — Les opéra- 
teurs d’un champ spinoriel quantifié peuvent se 
développer en ondes planes : 


V(r, £, m) — 2 f dfpiatt 9(p) ut+)(9(p, m) eilp.r—Et) 
$ 
+ a+ (p) (pp, m)eilp.r+Æ0)], (65a) 
Pr, t,— m) — 2 f'apiat-0(p) u(H)(s\(p, — m) ei(p.r—EÆt) 
6 


+ a (p) ut (sp, — m) eilp.r+Æ)], (65b) 
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Les a® (p) sont opérateurs d'absorption d’une 
particule et celui d’une antiparticule, et les 
a®®t (p) opérateurs d'émission correspondants. 
Ces opérateurs obéissent aux règles d’anticommu- 
tation 


{a (p}, att)(s)(p') = {a 9+(p), a (9+(p) ) 
{a (p), a+ (p’) } — 5(s, s’) Hp — p'). 


Nous examinons maintenant si W(r,{,m) et 
(r, lt, —m) peuvent être, à un facteur de phase 
près, invariants sous certaines opérations d’inver- 
sion. À Ce sujet, nous nous intéressons surtout aux 


inversions unitaires À ou AM, car les fonctions 


propres & peuvent être normées de façon uniforme 
par rapport à ces inversions comme nous l’avons 
montré aux paragraphes 3 et 4. 

Sous l’opération des inversions unitaires, nous 
supposerons que les opérateurs d’émission et d’ab- 
sorption du champ ne changent pas de phase, c’est- 
à-dire que ces opérateurs soient invariants d'image 
par rapport à ces inversions. Alors nous devons 
mettre (Cf. [1], $ 9) 


(ja (p) — (Ma 6)(p) = as (p), 
(Pja#®(p) = (PMja@(p) = a p),| 
(CTia& S(p) = (CTM)a HO (p) = amot(p), (109) 
(CPT)a)(s) (p) - Re Pte — as p; 
Cela étant, les opérateurs de champ W(r,t, m) 


et V(r,t,— m) seront invariants sous les inver- 
sions unitaires À et AM, si les conditions corres- 
pondantes sur les facteurs de phase € des fonctions 
propres 4 sont bien remplies (voir $ 3 et $4). Ainsi, 
pour toutes les inversions unitaires, nous avons 


CRETE) = (+ ou —) Wir,t,m), (674) 
CORNE mn EE our Em); (676) 
(AM) Wir, t, m) — (+ ou —) V(r,t, m), (67c) 
(CAM) (r,t, — 7m) = (H\ou —) Fir,t, — m) (674) 


Quelques cas particuliers de ces résultats, dont 
nous nous sommes servis dans des études anté- 
rieures [5], [6], trouvent ainsi leur démonstration. 

Nous nous intéressons surtout aux inversions 
CPT, M, CT, PM. Posons maintenant : 


Her, F m) — a) (s)(p) u(Æ)(s)(p, mn) eùtp.r+At) 


(684) 
(68b) 


+- a) (st (— P, mn) u(F)(s) (— D, In) eil-p.r+ÆEt), 


MS, t,m) — as) (p) uÆ(s) (p, m) eilp.rF£t, 
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Or, 4, m) = a (p) u HO (p, m) eilp.rF Et) 
+ aP6)+(p, m) uP(s(p, m) eilp.r£Et), (68c) 

LACNCS t,m) — as) {p) us) (p, m) ei(p.rF At) 
+ a (— p, m) uD OS (= p, m) eil-p.rFÆ), (68d) 


« Si nous choisissons des facteurs de phase € qui 


remplissent à la fois les conditions (63) et (64), alors, 
sous l’hypothèse (66), on montre que 


(CPT) FT LM) = EoUu—) Fr, €, m), (69a) 
(M) Fr, &, m) = (+ ou —) Fr, 1, m), (69b) 
(ÉTDIPO (nr, mon = Ir 5m), (69c) 

(PM) RACE lt, m) = (+ ou —) LAN t, Im). (694) 


Donc, au lieu de W(r, 4, m), mêmes les combi- 
naisons (68) sont déjà invariantes par rapport à 
CPT, M, CT, PM respectivement. 

On remarque finalement qu'aucune des inver- 
sions unitaires n’engendre le changement d’héli- 
cité, tandis que toutes les inversions antiunitaires 
engendrent un tel changement. 


9. Conclusion. — Par rapport à toutes les inver- 
sions unitaires, la normalisation uniforme des 
ondes planes monochromatiques est réalisable sous : 
deux possibilités différentes d’un signe seulement. 
Une normalisation uniforme est possible même 
pour plusieurs opérations d’inversion simultané- 
ment, en particulier pour CPT, M, CT, PM. 

Dans le cas des inversions antiunitaires, une 
normalisation uniforme est possible seulement 
quand la matrice de transformation correspon- 
dante Q est symétrique (C, PTM). Cependant, 
pour les inversions TM, CM, CPM,il peut y avoir 
une normalisation quasi-uniforme, normalisation 
qui est uniforme pour les états propres avec + m 
et aussi, mais différemment, pour les états propres 
avec — m. Pour T, PT, CP, la normalisation ne 
peut être uniforme mais, dans la meilleure possi- 
bilité, elle n’est différente que pour les états d’héli- 
cité différente. 

Grâce à l’existence de normalisation uniforme 
des ondes planes monochromatiques dans le cas 
des inversions unitaires, il existe des opérateurs de 
champ qui sont, au plus à un signe près, invariants 
par rapport à ces inversions. Les grandeurs (65) 
et (68) sont de tels opérateurs. 
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APPENDICE 


En adoptant une certaine représentation parti- remplissent les conditions (2) et (3). Prenons par 
culière des matrices de Dirac, nous donnons 1ci les exemple 
expressions explicites des fonctions propres 4, qui 


0 — 10 — I 0 0 À Fe 
Æ 4 Fe — -1 
où 
0 1 0 er 1 0 
— = — > A-2 
O1 [, de LE {; 0 À C3 L EL il ( ) 
Nous pouvons alors choisir pour les fonctions w, les spineurs suivants. 
— P(p + P3) 
(H)CHD LCD he — P(P1 + iP2) : 
uÿ (P, m) = uw (—p, — m) = N Fe RE : (A-3a) 
(E + m) (p1 + 1p2) 
— P(P1 — Po) 
(+1) 1) Re p(p + pa) à 
uo (p, m) — Uÿ fe P, m) = N Der (E a m) (Dre ipa) ; ‘ (A 3b) 
LT m)(p EPS 
[ _(£ + m){(p +p:) 
ug FD (p, HO) == DATE P, — m) = N rss , (A-3c) 
P\p + P3) 
Es P(p1 + ip2) 
UD —(£ + om) (p, — A 
PP er Dore Eee ne 2 . Pal (A-34) 
- P(Pp1 — : 
LEA SPP je, P3) 
P(p + Ps) 
HD p, m) = HD P(P1 + ip2) 
(72 ) = |) — N 
Ô Ù (P, — m) BR Aloe (A-3e) 
D (E + m){(p: Fip:) 
BE P(P1 — ip2) 
(HD, Re er : = Pipe 
uÿ ( DU — — À 3 
1 ) 0 (p In) PT Se . (pa y. ip») (A-3/) 
=, (Fm) (pp) 
TT (Em) = + pa) 
(EE) ; . CHA) (HE In) ENT D 
u D, m\ = es Je Pi P2) 2 
0 (— p, m) uQ (P, mm) “ue : te (A-3g) 
Le — p(p1 + 1P) 
Em) {pi —"ip:) 
C1) __, (1) : + m) (p + p 
up (— p, m) = u En) AN (Z ee 
0 (p mm) D ue ia (A-3h) 
! P(p + ps) 


Dans ces expressions, 
1 
N =; [pE(p + ps) (E + m)j—22. 


Manuscrit reçu le 11 février 1960, 
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LA DÉTERMINATION DE LA DIRECTION D'ARRIVÉE 
ET DE LA POLARISATION DES ATMOSPHÉRIQUES SIFFLEURS 


Deuxième partie : Les résultats. 


Par JEAN DELLOUE, 


Laboratoire de Physique de l'École Normale Supérieure, 


l Résumé. - Dans un article précédent on a décrit Pappareillage utilisé pour déterminer les direc- 

tions d'arrivée des atmosphériques siffleurs et leur polarisation. 
L'interprétation, assez délicate, des diagrammes observés sur les oscillographes de mesure fait 
l’objet du présent article. Une représentation simple des directions d’arrivée est proposée et les 
règles à observer pour obtenir des diagrammes de phase lisibles sont exposées. 

L'étude de la polarisation montre que, moyennant certaines hypothèses simplificatrices, les 
conditions moyennes le long de la trajectoire d’un siffleur peuvent être représentées par un seul 
paramètre + qui, avec l’angle entre la direction d’arrivée et le champ magnétique, détermine la 
polarisation dans le plan d’onde. La polarisation observée dépend alors de la direction d’arrivée 
et se déduit de la précédente par une projection cylindrique. Des diagrammes de polarisation cons- 
tituent la seule façon commode de remonter de la polarisation observée à celle dans le front d’onde 
et au paramètre +. 

L'exploitation du matériel représente une somme d’efforts considérables pour un pourcentage 

de succès faible. On a pu cependant mettre en évidence que la direction d'arrivée varie au cours 
d’un même siffleur en faisant un angle toujours assez faible, compris entre 5 et 259, avec la direc- 
tion du champ magnétique. Il n’a pas été mesuré d’angles inférieurs à 5°. Enfin, la direction d’arri- 
vée tourue autour du champ terrestre. 

La polarisation observée dépend à la fois de la direction d’arrivée et des condilions rencontrées 
le long de la trajectoire. Les résultats sont compatibles avec les directions d’arrivée déterminées 
mais conduisent à des densités électroniques plus faibles que ne lavait prévu Storey. Il ne semble 
pas dans ces conditions que l’approximation utilisée par Storey pour déduire les propriétés des 
sifileurs à partir de la théorie quasi-longitudinale de Booker soit suffisante. Un calcul plus complet 
rendrait peut-être compte de l'absence de directions de propagation très voisines de celle du champ 
terrestre. 

Bnfin, l’observation de suites brèves de quelques forts siffleurs présentant à l’intérieur d’une 
suite donnée les mêmes caractéristiques de polarisation et de déphasages conduit à supposer l’appa- 
rition subite et très localisée dans la très haute atmosphère de centres ionisés dont la durée de vie 
est faible, de l’ordre de quelques minutes. 


Abstract. — In a previous article an equipment for the determination of the direction of arrival 
and of the polarization of whistling atmospherics was described. 

The rather difficult interpretation of the patterns recorded is now discussed. A simple way of 
representing the directions of arrival is proposed and the rules which are to be observed to obtain 
usable phase patterns are explained. | non 

A study of the polarization of whistlers shows (hat, provided à few simplifying assumptions 
are made, the average conditions along the trajectory of a whistler can be represented by a single 
parameter + which, along with the angle between the direction of arrival and the magnetic field, 
determines the polarization ellipse in the wave front. The polarization observed depends then on 
the direction of arrival and can be deduced from the previous one by a cylindrical projection. 
Polarization diagrams are the only convenient way of solving the inverse problem and to determine 
the parameter +. és 

Experimentation on the field requires a considerable amount of effort and the percentage of 
success is low. It has been nevertheless possible to show that the direction of arrival rotates 

; around the magnetic field and makes a rather small angle, between 5 and 25°, with that direction. 
No angle smaller than 5° was observed. The projected polarization depends on the direction of 
arrival as well as on the conditions encountered along the propagation path. Our results are con- 
sistent with the directions of arrival observed, and lead to electronic densities smaller than those 
quoted by Storey. The approximation used by Storey Lo deduce the properties of whistlers 
from Booker’s quasi-longitudinal theory appears not to be sufficient. A more complete calcula- 
tion might account for the absence of direction of propagation in the vicinity of the magnetic 
field. ES 

The occurrence of sequences of several strong whistlers with the same polarization and phase 
characteristics within the same sequence could be explained in terms of the sudden appearance 
in the outer atmosphere of ionized centres with the short life time of about a few minutes. 


5. L'interprétation des mesures. 


5-1. Détermination des directions d’arrivée. — 
5-1.1. ASPECT GÉOMÉTRIQUE. — Considérons trois 
stations réceptrices situées aux sommets d’un 
triangle rectangle et mesurons les déphasages ©, 
®’ entre les signaux reçus d’une part à la station 
centrale O, au sommet du triangle, et d'autre part 
aux deux stations auxiliaires À, A’ (fig. 13). Ces 


Encre 


déphasages sont reliés aux angles 0, © qui définis- 
sent la direction d’arrivée de l’onde plane par 

D, = 2rd/JÀ d — OA 
D — 2nd'J\  d'— OA! 


D — ®, cos 0 cos 


L D’ — D; sin 0 cos w 


Nous supposons 0 < 6 < x/2, 0 pouvant prendre 
toute valeur. 
Ce syst ème peut encore s’écrire 


k == cos 0 cos © _ RAD] 
k' = sin 0 cos y RARDA ID 


RU EN: 


EN ME CES | 


D'où, en inversant 


tg 0 
cos © 


k'[k 
(GO 0) ETES 


Ces formules montrent que : 

a) si les déphasages sont mesurés expérimen- 
talement avec une erreur absolue A® donnée 
(comme c’est le cas), il y a intérêt à disposer de 
bases aussi longues que possible pour diminuer 
erreur relative AD/®. Cette longueur sera cepen- 
dant limitée par les caractéristiques des relais 
hertziens utilisés et le profil des liaisons, et par la 
nécessité de pouvoir considérer le front de l’onde 
comme homogène aux trois stations ; 

b) les erreurs sur 0 et © sont minima, pour une 
erreur absolue donnée sur les mesures de phase, 
lorsque 0 est nul. Si l’on a des raisons de considérer 
une direction d'arrivée comme probable, il y a 
intérêt à orienter l’une des bases au plus près du 
plan vertical contenant cette direction. Dans nos 


I 


I 
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mesures, l’une des bases est toujours le plus près 
possible de plan du méridien magnétique moyen. 


5-1.2. LA REPRÉSENTATION DES DIRECTIONS 
D'ARRIVÉE. — Nous avons adopté pour représenter 


les directions d’arrivée des siffleurs, un diagramme | 
en k, k! qui jouit de la propriété de donner immé- 


diatement l’angle 0 (formules IT), le point figu- 
ratif m n'étant autre que la projection sur le plan 
horizontal de M, extrémité du vecteur OM de 
module unité porté par la direction d’arrivée 
(fig. 13). Sur ce diagramme les lieux d’égale élé- 
vation © sont des cercles de centre O et les points 
figuratifs des directions faisant un angle constant 
avec une direction donnée sont situés sur des 
elbpses (fig. 14). En effet, supposons cette direction 
de référence OH dans le plan zOA, d’élévation D. 
L’angle (OM, OH) = e est donné par 


cos £ — cos 0 cos © cos D + sin 6 sin D 
d’où, d’après IT 
(k cos D — cos €)? — sin? D(1 — k? — k’?) 
soit 
(k — cos D cos €)? + k’? sin? D = sin? D sin? &. 


C’est l’équation d’une ellipse dont le centre a 
pour coordonnées (cos D cos e, 0) et dont les axes, 
parallèles à OA et OA’, ont pour mesures respec- 


ak 
fu k” 
D) y 

1 


FiG. 14, — Diagramme kk’ et ellipses à e constant. 


tives sin D sin & et sin e. Lorsque la direction de 
référence OH n’est pas dans le plan vertical zOA, on 
obtient je réseau d’ellipses correspondant à partir 
du précédent par une rotation d'angle © (@ est 
l’angle que fait la base Nord-Sud avec le méridien 
magnétique, lorsque la direction de référence OH 
est celle du champ magnétique local). Grâce à ce 
réseau d’ellipses « d’égal e », on peut rapidement 
apprécier la façon dont un siffleur arrive au voisi- 
nage du champ magnétique. 


EN°:7 ï 
9-1.3. L’APPARENCE DES SIFFLEURS SUR LES CA- 
THODIQUES DE PHASES. — Nous avons expliqué plus 
haut qu’un siffleur faisait apparaitre sur chaque 
cathodique de phase des segments (ou des ellipses 
très aplaties) dont l’angle avec la direction d’une 
des paires des plaques du cathodique donnait la 
mesure du demi-déphasage ®/2 correspondant. 
Cette paire de plaques est choisie de façon que ces 
segments soient lisibles, c’est-à-dire forment un 
angle aussi grand que possible avec l’axe de dérou- 
lement du film. Une appréciation sommaire des 
angles à attendre est donc nécessaire avant toute 
chose. Le plus grand nombre de mesures a été 
effectué dans la région parisienne avec les longueurs 
de bases suivantes : 


Nord-Sud 17,5 km soit D, — 1150 à 5 500 c/s 
Est-Ouest 22,5 km soit D, — 1500 à 5 500 c/s 


et un angle de 80 entre le base Nord-Sud et le méri- 
dien magnétique de la station centrale, la station 
Sud étant situé à l’Est de ce méridien. Pour la base 
Est-Ouest, la station centrale est à l’Est et la sta- 
tion auxiliaire à ’Ouest. Nous avons tracé ci-contre 
(fig. 15), pour le cathodique de phase Nord-Sud, les 


Phase 0 
appareillage 


Phase 0 
signaux 
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Fire. 15. — Aspects des cathodiques NS et EW 
et dispositions relatives des stations. 


directions correspondant à diverses élévations pour 
des directions d’arrivée situées dans le plan vertical 
de la base Nord-Sud (8 —0). On remarquera que la 
direction correspondant au déphasage 0 a été 
choisie horizontale pour des raisons de lisibilité 
mentionnées plus haut. La direction marquée 
® — 900 correspond à des signaux arrivant en 
phase aux deux stations, mais à un retard du 
signal de la station auxiliaire par rapport au signal 
central] de ® — 1150: eile fait donc l’angle 
D,/2 — 570 1/2, porté dans le sens des aiguilles 
d’une montre (sens des retards auxiliaires/central) 
à partir dela direction de déphasage nul. De même 
nous avons tracé pour le cathodique de phase Est- 
Ouest, les directions correspondant à divers azi- 
muths pour une élévation constante © = 650 (celle 
du champ magnétique terrestre). Cette fois, la lisi- 
bilité impose une direction de déphasage nul verti- 
cale et la direction marquée 0 — 09 fait un angle 
®,/2 — 750 avec la verticale. 


5-2, La mesure ‘et l’interprétation de la polari- 
sation. — 5-2.1. LA POLARISATION PRÉVUE PAR LA 
THÉORIE MAGNÉTO-IONIQUE. — La théorie magnéto- 
jonique donne le rapport À des composantes du 
champ magnétique d’une onde plane : 

{ |] située dans le plan contenant la normale au 
plan d’onde et la direction du champ ma- 
gnétique (plan magnétique). 

ki perpendiculaire à ce plan, 


R= nl = ir + 
y 1 
2 1—% 


.sin ec tg € 
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— angle de la normale au plan d’onde avec 
le champ magnétique. 

y — pufp — Hefpme rapport de la gyrofré- 
quence à la fréquence de l’onde. 

x = péfp? —4r Ne?]mp? carré du rapport de 
la fréquence critique à la fréquence de 
l'onde. 

Cette formule est valable Jorsqu’on néglige les 
collisions, approximation certainement correcte aux 
très grandes hauteurs où se localise la plus grande 
partie du trajet parcouru par les siffleurs. Lorsque 
x < 1 le rayon extraordinaire correspond, au 
signe + dans la parenthèse, tandis que si x > {il 
correspond au signe —. Le rayon ordinaire corres- 
pond à l’autre signe. Enfin, pour des directions de 
propagation faisant des angles aigus avec la direc- 
tion du champ magnétique, le rayon extraordinaire 
est polarisé dans le sens des aiguilles d’une montre, 
le rayon ordinaire dans le sens inverse. 

MO étant la direction d’arrivée (, 0) d’un sif- 
fleur et OH, située dans le plan zOA, celle du 
champ magnétique terrestre, deux polarisations 
sont possibles, représentées par deux ellipses dont 
les axes <ont situés l’un dans le plan magné- 


où 


[0] 


F1G. 16. — Polarisation dans le plan d'onde. 


tique MOH et l’autre dans le plan perpendiculaire 


(fig. 16). 


5-2.2. APPLICATION AU CAS DES SIFFLEURS. — 
La quantité y — pa/p — Hefpme dépend de la 
fréquence qui nous intéresse (5 500 c/s), et de la 
valeur du champ # sur la trajectoire du siffleur. 
Nous avons représenté sur la figure 17 l’allure d’une 
ligne de force du champ magnétique terrestre 
correspondant à une inclinaison de 659 à la surface 
de la terre, Le calcul présenté en appendice donne 
la construction de cette ligne et le tableau I les 
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valeurs des diverses grandeurs pour quelques 
altitudes. 


TABLEAU T 


ALTITUDE (km) (] 4700 3900 5700 7 000 
Champ (gauss) 0,46 0,2 0,1 0,05 0,03 
Gyrofréquence ” 
(kc/s) 1,290 560 280 140 84 
y 235 102 51 26 16 


La quantité x — p$/p? = 4x Ne?[mp? (4 2,5 N 
à 5 500 c/s) peut varier de 0 à 105 suivant les den- 
sités rencontrées le long de la trajectoire. Comme 
les dispersions observées obligent à admettre, dans 
le cadre de la théorie de Storey, des densités élec- 


F1G. 17. — Ligne de face du champ terrestre 
correspondant à la région parisienne. 


troniques assez élevées sur la plus grande partie 
de la trajectoire, nous pouvons admettre à coup 
sûr que x > 1 (N > 1/3): en fait, si les valeurs 
de N avancées par Storey sont réelles (W — 400), 
æ est de l’ordre de 1 200. L'expression précédem- 
ment écrite pour f peut alors se simplifier : 


Î#— (Tr{2) sin e tg « avec + = yJx 


et f est négatif tant que la direction de propagation 
fait un angle aigu avec celle du champ magnétique. 

La quantité T, que nous venons d'introduire, peut 
aussi S’écrire : 


5 = ylx = (HIN).(f{2ec) = 180H/N à 5 500 cjs. 


On peut donc, dès maintenant, faire les re- 
marques suivantes : y est inférieur à 200 sur la 
plus grande partie de là trajectoire ; x est supérieur 
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à 100 dès que N = 75 électron/em$. D’autre part, 
si e << 500, 


sin etg € < 1. 


Donc, dans cés conditions, f < 1 et R > 0,5. 
Il faut donc s’attendre à une polarisation fortement 
elliptique. Si x est dix fois plus grand, la polari- 
sation est pratiquementtoujours circulaire, quel que 
soit <. Mais ceci ne signifie d’ailleurs pas que la 
polarisation en projection observée soit nécessai- 
rement aussi forte, car elle dépend de la direction 
d'arrivée. 


5-2.3. LA GRANDEUR DE À DANS LE CAS DES 


SIFFLEURS. — Il suffit de calculer les valeurs de R 
comprises entre 0 et 1, c’est-à-dire d’évaluer 
VOTE | à en ca D  E 
Posante— cotg ©, (0 © <"x)2) 1 vient 


ft {p.0 7/2. 
Nous tracons les courbes ÆA(<e) pour diverses 
valeurs de + par le processus de caleul 


Te. re) 
e—> sin e g € — cotg © — ow * 18 = R 


De l’examen de ces courbes (fig. 18) on peut 


conclure que pour des directions faisant moins 
de 600 (par exemple) avec la direction du champ 
magnétique terrestre, la polarisation est prati- 
quement toujours circulaire dès que tee 0, 
Ce résultat avait été prévu dès la fin du para- 
graphe précédent. 


5-2.4. LA POLARISATION DES SIFFLEURS ET LE 
raPport JAN. — I] faut remarquer que la pola- 
risation d’un siffleur dépend des valeurs de la den- 
sité électronique et du champ magnétique tout le 
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long de la trajectoire. Il est-pratiquement impos- 
sible de calculer exactement l’effet d’un tel trajet, 
ne serait-ce d’ailleurs que parce que la loi de varia- 
tion de nous est inconnue. Si l’on admet comme 
hypothèse de travail que la quantité H/JN est sen- 
siblement constante le long de la trajectoire (car il 
est vraisemblable que N décroît, comme Æ, quand 
l'altitude croît), on peut calculer cette quantité à 
partir de la polarisation observée. Si, par exemple, 
H = 0,05 et N — 450, on trouve pour + la 
valeur 0,02 : La polarisation résultante serait donc 
pratiquement circulaire quel que soit l’angle e. 


5-2.5. POLARISATION DANS LE PLAN D'ONDE ET 
POLARISATION EN PROJECTION. — L'emploi de 
cadres verticaux pour les mesures de polarisation 
revient à ne considérer que les composantes hori- 
zontales du champ magnétique de l’onde. Il est 
d’ailleurs à remarquer que, le sol pouvant être 
considéré comme parfaitement conducteur aux fré- 
quences très basses étudiées, il n’y a pas de compo- 
santes verticales décelables, et les composantes 
horizontales sont doubles de celles de l’onde inci- 
dente loin du sol. Autrement dit, l’ellipse observée 
sur le cathodique de polarisation est la projection 
sur le plan horizontal de l’ellipse représentant le 
champ magnétique dans le plan d’onde. On voit 
donc que pour une même polarisation dans le plan 
d’onde, et des directions d’arrivée différentes, l’as- 
pect de l’ellipse en projection sera très variable. 

Nous avons reporté ci-contre sur des dia- 
grammes Xk', en les centrant sur le point m repré- 
sentatif, des ellipses semblables à celles qu’on 
observerait pour les directions d’arrivée corres- 
pondantes (fig. 19). Chacun des réseaux est tracé 
pour une valeur donnée de 7 : les ellipses tracées en 
tirets correspondent à des ellipses du plan d’onde 
dont le grand axe est dans le plan magnétique 
(rayon extraordinaire), les ellipses en trait plein aux 
ellipses du plan d’onde dont le grand axe est per- 
pendiculaire au plan magnétique (rayon ordinaire). 

Ces diagrammes ont été obtenus par une méthode 
approchée au moyen d’une projection optique 
d’ellipses tracées sur plexiglass et convenablemeut 
orientées par rapport au plan de projection. 

On remarquera sur ces réseaux que l’orientation 
de l’ellipse projection et son excentricité dépendent 
de façon très complexe de l’excentricité de l’ellipse 
du plan d’onde et de la direction d’arrivée de cette 
onde. Mais il est clair qu’étant donnée une direction 
d'arrivée déterminée, une ellipse quelconque en 
projection ne correspond pas nécessairement à une 
ellipse « réelle » du plan d’onde associé : 1l y a une 
relation à respecter entre l’orientation et l’excen- 
tricité. La mesure de la polarisation apporte donc 
en plus d’un renseignement sur l’état de polari- 
sation de l’onde reçue (c’est-à-dire sur +) une con- 
firmation de l’exactitude des mesures de direction 
d’arrivée. 
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6. La mise en œuvre du matériel. 
Les résultats. 


6-1. La mise en œuvre du matériel. — Le dé- 
ploiement de matériel et de personnel nécessité par 
. les mesures est considérable. Chaque station auxi- 

haire exige au moins un technicien capable d’ériger 
les antennes, d'établir la liaison hertzienne et de 
s'assurer du fonctionnement correct de l’appareil- 
Jage. A la station centrale, outre les deux liaisons 
hertziennes à établir et à maintenir et la surveil- 
lance d’un ensemble de mesures très complexe, il 
faut assurer la photographie continue des oscil- 
lographes (charge et décharge de la caméra toutes 
les vingt minutes), le contrôle de la qualité des 
prises de vues (développement des films de 16 mm 
au fur et à mesure, nécessitant un laboratoire spé- 
clal parfaitement équipé, et une forte pratique du 
traitement des émulsions rapides), le repérage des 
siffleurs et l’enregistrement magnétophonique syn- 
chronisé sur la prise de vues. Un minimum de deux 
personnes s’impose. Ainsi chaque séance de mesures 
mobilise de quatre à sept personnes. 

Les siffleurs se manifestant principalement la 
nuit, ce personnel travaille dans des conditions 
particulièrement pénibles. D’un autre côté, il est 
pratiquement impossible de prévoir les nuits où des 
siffleurs se manifesteront, et avec une intensité 
suffisante pour que les mesures aient une valeur : 
les séances fructueuses sont donc rares (en moyenne 
une sur cinq). On cherche évidemment à éviter les 
déploiements inutiles en n’opérant que pendant 
les périodes d’activité (liées d’ailleurs à l’activité 
solaire), mais même dans ces périodes les trous sont 
nombreux. Enfin, les heures d’apparition des 
siffleurs sont très variables, sans qu’on puisse même 
être sûr de l’heure moyenne à attendre. Or, dans 
la plupart des cas, au moins une des stations auxi- 
liaires fonctionne sur accumulateurs lui assurant 
une autonomie de 4 à 5 heures, après quoi une 
recharge d’une dizaine d'heures est nécessaire avant 
la séance suivante. Le facteur chance a donc une 
importance considérable. 

A ces difficultés imputables aux phénomènes 
observés s'ajoutent celles provenant des ennuis de 
matériel (heureusement rares), du temps (lhiver 
est la période la plus favorable) en particulier des 
problèmes causés par le vent dans les aériens qui 
sont obligatoirement en des emplacements tres 
dégagés (phares, falaises, toitures), des particu- 
larités de la propagation aux fréquences élevées 
utilisées (surtout gênantes lors de l’établissement 
des liaisons, quoiqu'il soit aussi arrivé de voir des 
liaisons déjà éprouvées se montrer subitement 
impraticables). Deux emplacements ont été utilisés 
avec desmodes d'exploitation totalement différents. 

— Les premiers essais ont été faits dans la baïe 
du Mont Saint-Michel, la station centrale se trou- 
vant à Granville, les stations auxiliaires à Carteret 
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et au Cap Fréhel. Les bases correspondantes sont 
très longues (60 km), et les trajets par la route entre 
la station centrale et les stations auxiliaires assez 
grands (Carterat 90 km, Fréhel 140 km) empê- 
chaient d'envisager à chaque séance un déploiement 
à partir de la station centrale. Nous avons donc 
fonctionné chaque nuit pendant une période de 
deux semaines, les opérateurs vivant dans le voisi- 
nage de leur station. D’autre part, l'éloignement de 
Paris rendait impensable toute modification un peu 
conséquente de l’appareillage. Lorsque à la fin des 
deux semaines des siffleurs se sont enfin manifestés, 
il est apparu, outre un certain nombre de résultats 
fondamentaux, la nécessité de revoir la conception 
de l’ensemble de mesures (voir 4 partie). Le maté- 
riel a été alors rapatrié dans la région parisienne où 
toutes les mesures ultérieures ont eu lieu. 

— Le second emplacement, quoique ayant 
l'avantage considérable d’être près de Paris, a un 
certain nombre de défauts :Jes bases sont beaucoup 
plus courtes (17,5 et 22,5 km), les stations auxi- 
liaires doivent être remises en place à chaque 
séance et démontées à la fin, le voisinage des lignes 
de distribution force et de nombreux appareils 
électriques causent un brouillage des fréquences 
basses et des parasites impulsifs. C’est surtout 
l'aspect précaire des stations auxiliaires qui est le 
plus gênant : en effet, à partir du moment où une 
séance est décidée, il faut compter environ une 
heure (trajet et mise en place) avant le début des 
mesures. Or, ces décisions sont souvent prises 


- lorsqu’on entend les premiers siffleurs et qu’on juge 


leur nombre et leur force suffisants. Plus d’une fois, 
les siffleurs avaient disparu lorsque les mesures ont 
pu commencer. L’une des deux stations auxiliaires 
est équipée d’une camionette, principalement pour 
gagner, grâce à une plate-forme de toit, les quelques 
décimètres nécessaires pour assurer la liaison cor- 
recte Est-Ouest. L’autre station est transportée 
dans une voiture de tourisme. 


6-2. L'enregistrement synchrone sur magné- 
tophone. — En même temps que la caméra photo- 
eraphie les écrans des oscillographes cathodiques, 
un magnétophone enregistre les atmosphériques à 
la sortie d’un des amplificateurs large bande de la 
partie polarisation. 

Pour permettre l’identification de chaque siffleur 
sur le film et sur la bande magnétique, l’opérateur 
qui observe les écrans des cathodiques ferme un 
contact chaque fois qu’un siffleur provoque une 
déflection notable. Ce contact allume d’une part un 
témoin néon placé entre les oscllographes et 
d'autre part provoque l’enregistrement d’un top 
à 400 périodes sur la demi-piste inférieure de la 
bande magnétique (la demi-piste supérieure étant 
réservée à l’enregistrement des phénomènes). 

Lorsqu'on examine le film développé on peut 
amsi prélever très rapidement les siffleurs. D'autre 
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part, à la lecture de la bande magnétique les tops 
400 périodes actionnent un voyant qui permet de 
sélectionner les parties utiles de l’enregistrement. 
Enfin, ce top sert de repère précis pour déterminer 
à quelle portion exacte du siffleur correspond 
chaque « tranche » photographiée sur les catho- 
diques, comme nous allons le préciser. 


6-3. L’analyse spectrale des siffleurs. — L’enre- 
gistrement magnétique des siffleurs Jugés intéres- 
sants est repris dans un analyseur basse fréquence 
à tambour magnétique (Sonagraph) qui trace auto- 
matiquement un diagramme portant le temps en 
abscisses, la fréquence en ordonnées (de 50 à 
8 000 c/s), un noircissement plus ou moins fort tra- 
duisant l’importance de chaque composante spec- 
trale à un instant déterminé. On obtient en une 
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dizaine de minutes une représentation très com- 
mode du phénomène qui permet de mesurer la dis- 
persion, de voir si l’on a affaire à un siffleur long ou 
à un siffleur court, et de dépister les siffleurs d’éclair 
multiple. 

Si par un couplage électrique, on injecte dans 
l'enregistrement magnétique le signal 400 périodes 
de repérage, celui-ci se traduit sur le diagramme 
spectral par un ou plusieurs tirets horizontaux. 
Une simple mesure de longueur permet de repérer 
alors sur le diagramme les « tranches » photo- 
graphiées : avantage considérable, par exemple 
quand la «largeur » du siffleur a été si faible qu’une 
seule tranche a été photographiée et qu’on hésite à 
l'identifier sur le film, ou bien encore lorsque des 
atmosphériques importants semblent avoir brouillé 
les oscillogrammes, car un simple coup d’œil permet . 
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Fi@. 20, — Film d’unsiffleur et spectrogramme associé. 


de repérer la barre verticale correspondant à un 
tel atmosphérique. Nous avons tenté d'illustrer ces 
remarques sur la figure 20, 


6-4. Résultats obtenus avec la première version 
de l’ensemble de mesures. — La première version a 
été essayée dans la baie du Mont Saint-Michel, sur 
des bases de 60 km. Comme nous l’avons expliqué 
dans la quatrième partie, cette version avait deux 
défauts : faible sensibilité et rotations de phases 
parasites introduites par les circuits des relais 
hertziens, Malgré le manque de sensibilité, un bon 
nombre de siffleurs a pu être observé. Les mesures 
de polarisat on, non perturbées dans cette version, 


étaient utilisables quoique peu précises. Quant au 
défaut de rotation parasite, nous avons dit qu’il 
faisait apparaître à la place d’un segment une tache 
de contour elliptique dont la direction de plus 
grande extension donnait une mesure correcte du 
déphasage : nous avons done pu utiliser quand 
même les mesures de déphasages, en nous conten- 
tant d’une précision moins bonne. Nous avons 
obtenu les résultats suivants : 


6-4.1. LE SIFFLEUR, PHÉNOMÈNE COMPLEXE. — 
Nous avons tout d’abord constaté qu’on ne pouvait 
pas parier de la direction d’arrivée d’un siffleur. En 
effet, au cours d’un même siffleur, de largeur natu- 


N°7 


_relle suffisamment grande pour qu’on observe plu- 


sieurs «tranches » successives sur le film, les carac- 
téristiques de phase et de polarisation varient 


considérablement. Les segments de déphasage 


Ctournent », indiquant les déphasages qui varient 
dans le temps, et par conséquent une direction 
d'arrivée non fixe. Nous avons évidemment été 
intrigués par ce phénomène et c’est en cherchant la 
cause dans l’appareillage que nous avons mis en 
évidence les rotations parasites dues aux relais 
hertziens : mais comme expliqué plus haut, ces 
rotations n’invalident pas complètement ces me- 
sures. D'ailleurs, si le moindre doute subsistait à 
cesujet, l’observation que nous avons faite delois de 


- variation des déphasages très variées et même 


parfois pour quelques siffleurs de déphasages cons- 
tants, suffirait à montrer que l’appareillage ne peut 
être seul mis en cause. 

D'autre part, les ellipses de pclarisation obser- 
vées, quoique inutilisables parce que trop petites 
pour des mesures plus précises, ont révélé une 


variété de formes et d’orientations incompatibles 


avec des directions d'arrivée fixes. Il ne saurait être 
question d’invoquer le facteur + pour expliquer 
cette variation, car on ne voit pas bien pourquoi à 
une direction d’arrivée fixe, et par suite à un 
trajet déterminé dans la haute atmosphère, com- 
prendrait un +aussirapidement variable. D’ailleurs 
les variations d’orientation, qui, comme nous 
Pavons dit dans la cinquième partie, sont très liées 
aux variations de la direction d’arrivée, confirment 
à elles seules les mesures de déphasages. 


6-4.2. LA LARGEUR NATURELLE DES SIFFLEURS. 
— Une explication intéressante de la «largeur » des 
siffleurs, est la multiplicité des return-strokes dans 
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les éclairs [8]. Cette interprétation est suggérée par 
l’observation dans les oscillogrammes des siffleurs 
de battements dûs à des interférences de fréquences 
variées présentes simultanément. Ces battements 
peuvent se justifier en admettant le chevauchement 
de siffleurs composants de mêmes caractéristiques 
qui se suivent à des intervalles de 30 à 50 milli- 
secondes, ordre de grandeur de l'intervalle entre 
les return-strokes d’un même éclair. 

Nos observations nous conduisent à proposer 
comme autre explication possible de la largeur des 
siffleurs les arrivées successives de siffleurs élémen- 
taires provenant d’une seule et même décharge, 
mais ayant parcouru des chemins voisins suffi- 
samment différents pour justifier cet étalement dans 
le temps. La vitesse moyenne de propagation des 
siffleurs étant de l’ordre de 15 000 km/s et la lon- 
gueur de leur trajectoire d'environ 25 000 km à nos 
latitudes, on voit que les intervalles de 30 à 50 milli- 
secondes précédemment mentionnés peuvent s’ex- 
pliquer à partir de différences de longueur de tra- 
jets de l’ordre de 500 km, c’est-à-dire de 2 à 3 % du 
parcours total. Nous montrons en appendice qu’un 
ordre de grandeur de l’allongement de la trajectoire 
lorsque celle-ci est légèrement différente de la ligne 
de force et qu’elle aboutit à une certaine distance d 
du pied de cette ligne de force, est de 750 km pour 
d — 100 km. Or nous avons vu que les observations 
de Storey, et d’autres plus récentes d’ailleurs [9], 
montrent qu’on peut capter le même siffleur dans 
un rayon de plusieurs centaines de kilomètres. 

Nous pensons que cette remarque.et les résultats 
de nos mesures de direction d’arrivée et de polari- 
sation suggèrent comme explication dela «largeur » 
naturelle des siffleurs l’existence d’un faisceau de 
trajectoires voisines empruntées par l'énergie 
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autour de la ligne de force moyenne et que c’est 
les différences de longueur de ces trajectoires qui 
sont cause de l’étalement du siffleur dans le temps. 


6-4.3. LA ROTATION DE LA DIRECTION D'ARRIVÉE 
AUTOUR DU CHAMP MAGNÉTIQUE. — S1 nous portons 


sur un diagramme k#k’ les points représentatifs des . 


tranches successives d’un siffleur, nous obtenons un 
lieu qui présente les particularités suivantes : 1l 
s'approche rarement à moins de 5° du point À figu- 
rant la direction du champ terrestre ; 1l s’écarte 
rarement à plus de 250 du point À ; le plus souvent, 
il semble s’enrouler plus ou moins grossièrement 
autour du point k, avec parfois rebroussement de 
chemin. 

Nous avons dessiné ci-contre quelques lieux 
observés par nous et qui illustrent assez bien ce qui 
précède (fig. 21). 

Ces résultats peuvent s’interpréter comme suit : 
un siffleur peut être considéré comme un ensemble 
de paquets d’ondes, ou « siffleurs élémentaires ». 
Chaque siffleur élémentaire voyage le long de sa 
trajectoire propre, tandis que les ondes planes qui 
le constituent se propagent dans des directions 
faisant des angles # faibles avec la trajectoire qui 
est pratiquement parallèle à la direction moyenne 
du champ magnétique terrestre. 

Si donc nous menons par la direction d'arrivée 
que nous avons déterminée pour une tranche de 
siffleur un plan P parallèle à la direction du champ 
terrestre D, la trajectoire du siffleur élémentaire est 
une parallèle à à cette dernière direction située 
quelque part dans ce plan (fig. 22). La direction 


5" D Ne 


d'arrivée des ondes fait avec cette trajectoire 
l’angle e que nous évaluons sur le diagramme kk'. 
D’après nos mesures, cet angle e est compris 
entre 59 et 259, La théorie magnéto-ionique des 
siffleurs est donc confirmée dans la mesure où 
Pangle e est faible, mais nous avons vu dans la 


première partie que les surfaces de vitesses radiales 


ne laissent pas prévoir de singularités pour les di- 
rections voisines de celle du champ magnétique. 


6-4.4. L'IONISATION DANS LA TRÈS HAUTE ATMO- 
SPHÈRE. — I] nous a été donné d’observer plusieur; 
fois, lors de nuits où les siffleurs notables étaient 
rares, mais où l’on percevait un fond continu de 
siffleurs extrêmement faibles, des suites brèves de 
quelques beaux siffleurs présentant à l’intérieur 
a’une suite donnée les mêmes caractéristiques de 
polarisation et de déphasages. 

Nous proposons comme explication de cet inté- 
ressant phénomène l’apparition subite et très loca- 
lisée de centres ionisés dans la très haute atmo- 
sphère. La présence de siffleurs très faibles montre 
en effet que les éclairs-source étaient nombreux, et 
que seule la propagation s’opposait à la formation 
de siffleurs notables. La faible étendue de la zone 
ionisée imposait un pinceau de trajectoires bien 
déterminées aux siffleurs éventuels, ce qui explique 
l’identité des caractéristiques et de leur variation 
dans le temps. La durée de vie de ces centres 1onisés 
est faible, de l’ordre de quelques minutes. 


6-5. Résultats obtenus avec la dernière version 
de l’ensemble de mesures. — La dernière version de 
l’ensemble a été utilisée dans la région parisienne 
sur des bases de l’ordre de 20 km. Grâce à une 
sensibilité grandement augmentée, elle a permis des 
mesures plus précises, surtout en ce qui concerne la 
polarisation. La suppression de l’effet de rotation 
parasite a rendu les mesures de déphasages complè- 
tement valables. Malheureusement, après la mise 
au point de cet ensemble, il n’y a eu que de très 
rares apparitions de siffleurs, et toujours de faible 
durée (une nuit). Nous n'avons donc eu que peu 
d’enregistrements. Les résultats suivants ont pu 
cependant être obtenus. 


6-5.1. DIRECTIONS D'ARRIVÉE. — [L’allure géné- 
rale de la variation de la direction d’arrivée dans le 
temps a été confirmée, ainsi que les faibles valeurs 
de l’angle d’arrivée. Il n’a toujours pas été observé 
d'angle inférieurs à 5 degrés. 


6-5.2. SIFFLEURS A ÉVOLUTION RAPIDE. — Des 
variations extrêmement rapides de la direction 
d'arrivée (confirmées d’ailleurs par des variations 
correspondantes de la polarisation en projection) 
ont été observées : une « tranche » de siffleurs se 
résolvant parfois en un écheveau d’ellipses toutes 
distinctes avec passage continu d’une ellipse à la 
suivante. La direction d’arrivée ne peut pas alors 
être considérée comme fixe pendant un temps aussi 
bref que le millième de seconde, Nous pensons que 
ce phénomène se présente pour des positions parti- 
culières des stations de mesures par rapport à la 
gerbe des trajectoires des siffleurs élémentaires, 
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L'interprétation des enregistrements est alors très 
complexe puisqu’il faut associer correctement deux 

à deux les divers déphasages observés, On y arrive 

en passant par continuité d’un motif à l’autre et en 
_ Saidant de la polarisation pour éliminer les orien- 
tations incompatibles. Les caractéristiques de l’en- 
registrement photographique (vitesse du film et 
modulation de brillance) étaient mal adaptées à ce 


genre de phénomène. 
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6-5.3. LA POLARISATION. — La polarisation 
observée est toujours fortement elliptique. L’orien- 
tation et l’excentricité des ellipses sont en général 
compatibles avec la direction d’arrivée par les 
mesures de déphasages, des cas douteux se pro- 
dusant parfois pour les variations rapides men- 
tionnées plus haut, 

Nous avons représenté ci-contre un magnifique 
siffleur tournant à évolution rapide (fig. 23) en 


chaque point duquel nous avons reporté les ellipses 
de polarisation correspondantes. On pourra voir par 
comparaison avec la planche d’ellipses de polari- 
sation + —3 que les formes de ces ellipses con- 
firment bien les mesures de direction d’arrivée : 
c’est en fait un cas typique où la polarisation a 
permis de lever des doutes sur les associations des 
déphasages deux à deux. 


6-5.4. LE + ÉQUIVALENT. — Le peu de résultats 
en notre possession indique nettement un + équi- 
valent plus élevé que ne le laissent prévoir les con- 
séquences de la théorie magnéto-1onique que nous 
avons développées dans la cinquième partie. Nous 
avons observé couramment des ellipses dont le 
rapport des axes est 1/3 ou moins, pour des angies € 
de l’ordre de 209, ce qui conduirait à un + de l’ordre 
de 3, c’et-à-dire à des densités électroniques de 
l’ordre de 30 électrons par cm, 10 fois plus faibles 
que celles proposées par Storey pour expliquer la 
dispersion des siffleurs. 

Pour évaluer cette dispersion, Storey admet 
comme nous un +t représentatif des conditions 


moyennes sur la trajectoire, mais l'expression qu’il 
obtient n’est valable que si + est nettement plus 
petit que 1, ce qui n’est pas le cas pour notre esti- 
mation de la polarisation. 

Cette limitation sur + est d’ailleurs aussi une 
condition de validité des résultats obtenus par 
Storey à partir de l’approximation quasi-longitu- 
dinale, et exposés à la fin de notre première partie. 
Il est donc possible que les calculs de Storey doivent 
être repris plus en détail sans la simplification T1. 
Mais ce travail sort du cadre que nous nous sommes 
fixé. 


6-5.5. LES VARIATIONS DE T AU COURS D'UN 
MÊME SIFFLEUR. — Enfin les mesures de polari- 
sation ont montré que le + équivalent pouvait 
varier non seulement d’un siffleur au suivant, mais 
encore d’un siffleur élémentaire à l’autre au cours 
d’un même siffleur. Ce dernier résultat, rapproché 
des remarques que nous avons faites au sujet des 
ionisations fugitives et très localisées de la très 
haute atmosphère, semble indiquer que des inhomo- 
généités très grandes de la densité électronique se 
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produisent dans ces régions, au point que des tra- 
jectoires relativement voisines rencontrent des con- 
ditions très différentes. 


7. Conclusion. 


Nous nous étions assignés deux objectifs au com- 
mencement de cette étude : déterminer l’angle 
d’arrivée des siffleurs, et vérifier si la polarisation 
mesurée est compatible avec une propagation dans 
un milieu 1onisé en présence d’un champ magné- 
tique. 

Les observations ont très rapidement montré que 
le problème est complexe. On ne peut parler de la 
direction d'arrivée d’un sifileur, ce dernier étant 
constitué par une succession de siffleurs élémen- 
taires arrivant en apparence de directions diffé- 
rentes. Autre particularités inattendues : si les 
angles de direction d’arrivée avec le champ magné- 
tique sont toujours assez faibles, on n’observe pas 
d’angles inférieurs à 5° et ces directions tournent 
autour du champ magnétique terrestre. Nous sug- 
gérons comme explication de ces phénomènes l’exis- 
tence d’un faisceau de trajectoires voisines em- 
pruntées par l’énergie autour de la ligne de force 
moyenne : les différences de longueur de ces tra- 
jectoires sont cause de l’étalement du siffleur dans 
le temps et leur disposition spatiale par rapport à 
l’ensemble des stations des mesures est responsable 
de la variation de la direction d’arrivée. L'énergie 
se propage par rapport à ces trajectoires, qui sont 
d’ailleurs sensiblement parallèles au champ ter- 
restre local, dans des directions faisant avec celle-ci 
un angle toujours faible compris entre 59 et 250, 

Les mesures de polarisation confirment dans 
l’ensemble les déterminations de directions d’arri- 
vée, mais l’observation fréquente de rapports des 
axes des ellipses de polarisation inférieurs à 1/3 pour 
des angles plus petits que 259 semble indiquer des 
densités électroniques dans la très haute atmo- 
sphère environ 10 fois plus faibles que celles pro- 
posées par Storey pour expliquer la dispersion des 
siffleurs. 

Il paraît donc nécessaire de reprendre le déve- 
loppement de lapproximation quasi-longitudinale 
de Booker dans le cas des siffleurs avec des hypo- 
thèses plus générales que celles de Storey (x/y > 1), 
afin d'essayer de justifier les dispersions observées 
à partir de densités électroniques plus faibles. Il est 
possible qu’on mette alors en évidence des singu- 
larilés pour les directions de propagation très 
proches de celle du champ terrestre, comme parais- 
sent l’indiquer nos mesures. 

Nous n'avons pas pu, par suite des circonstances 
défavorables, adjoindre à ce travail une partie sta- 
tistique basée sur un très grand nombre d’obser- 
vations. C’est un de nos objectifs futurs, ainsi que 
la modification de nos procédés d’enregistrement 


pour les rendre utilisables avec les siffleurs à varia- 
tion rapide. 
Enfin, des résultats intéressants ont été obtenus 
sur l’aspect de l’ionisation dans la très haute atmo- 
sphère, mettant en évidence l’apparition de zones 
d’ionisation fugitive de très faible extension. 


APPENDICE 
Sur le champ magnétique terrestre 


La distribution du champ à la surface de la 
terre est approximativement celle d’un aimant 
placé au centre de la terre. Les composantes radiale 
et tangentielle du champ étant 


He 2 Micos 0) Hy = M sin 0/r 
l’angle du champ magnétique avec le rayon vecteur 
est donné par 


0 EG 
Res se A 


Par intégration on obtient l’équation polare 
d’une ligne de force : 


PI=NT) SUN A0 


r, Correspondant à l’apex de cette ligne de force. 
L’inclinaison étant de 659 à Paris on calcule 
ro = 2,15 À, R étant le rayon terrestre. La ligne de 
force correspondante a été tracée (fig. 17) et on y a 
porté la valeur de /{ en divers points, suivant la 
Loi : 
H = (Mr) (3 cos? 0 + 4)/2, 


On calcule aisément la longueur de la portion de 
ligne de force située au-dessus d’une altitude 
donnée r : c’est 


S = (ry[V3) (2 E : sh 22) 
avec 


shæ — V3 cos 0 et r = r, sin? 6. 


La longueur totale du parcours d’un siffleur à la 
latitude de Paris s’obtiendra pour r —= AR 
ro — 2,15 À, et elle est voisine de 25 000 km. 

On remarque que Æ varie relativement peu 
autour de 0,05 gauss au-dessus de r = 1,6 R : la 
longueur du parcours effectué dans ces conditions 
est de 15 000 km, soit les 3/5 du parcours total. 
L'hypothèse simplificatrice de la cinquième partie 
s'appuie sur cette considération. 

Sans chercher à calculer la longueur exacte des 
trajectoires possibles des siffleurs, on peut estimer 
grossièrement la longueur d’une telle trajectoire 
partant du pied austral de la ligne de force de réfé- 
rence et arrivant à une colatitude 0 dans l’hésmi- 
sphère boréal. Cette longueur / est comprise entre 


ë 
N°7, 


Ja longueur ZL de la ligne de force de référence et 
celle L' de la ligne de force arrivant au même 
point que cette trajectoire voisine : nous pren- 
drons ! —(L + L')/2. Ces deux lignes de force 
étant très voisines, on obtient leur différence de 
longueur en différentiant l’expression de l’are déjà 
obtenu. Il vient 


AS = — 2 (sin 0)—1{($S cos 8 + R ch x) d8. 


Ce qui donne, à la latitude de Paris, dS — 
1 500 km/degré, soit pour la trajectoire du siffleur 
dS — 750 km/degré. 


Autrement dit, pour deux siffleurs élémentaires, 
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issus d’un même éclair mais arrivant par deux tra- 
Jectoires voisines en deux points du même méridien 
distants d'environ 100 km, il y a une différence de 
parcours d'environ 750 km qui correspond pour 
une vitesse de propagation de 15 000 km/s à un 
décalage dans les temps d'arrivée de 50 milli- 
secondes, Le processus d’élargissement suggéré 
dans la sixième partie permet donc d’expliquer la 
longueur observée pour les siffleurs en postulant 
des zones d’arrivée ayant des dimensions de l’ordre 
de quelques centaines de kilomètres, ce qui est 
conforme aux observations de Storey et autres [9]. 


Manuscrit reçu le 20 novembre 1959. 
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CALCUL DES SUSCEPTIBILITÉS MAGNÉTIQUES PRINCIPALES 
DE LA MOLÉCULE DE BENZÈNE 


Par Jean BAUDET, JEAN GUY et Jacques TILLIEU, 


Laboratoire de Chimie Physique de la Faculté des Sciences de Bordeaux 
et Groupe de Chimie Théorique de l’École Polytechnique. 


Résumé. La méthode d’approximation, que nous avions précédemment utilisée pour le 
calcul des susceptibilités diamagnétiques de nombreuses molécules organiques On aromatiques, 
est appliquée au benzène. Tandis que les valeurs théoriques et mesurées sont proches (écart relatif 
de 7 %) pour un champ dirigé perpendiculairement au plan de la molécule, on corstale par contre 
une différence plus sensible pour lPautre susceptibilité principale (34 % environ). Un calcul plus 
élaboré permet d’amener y, à 6 °/ de expérience mais ne conduit à aucune amélioration pour yn- 
Malgré tout, ces résultats apparaissent comme les meilleurs obtenus jusqu’à présent à notre 


connaissance. 


Abstract. — The approximation method we previously used for the calculation of the diamagne- 


tic susceptibilities of many non-aromatic organic molecules, is now applied to benzene. 


The theore- 


tical and experimental values are quite close (relative difference of 7 %) when the field direction 


is normal to the molecular plane ; 


on the other hand, a greater discrepancy is observed for the 


other principal susceptibility (about 34 %). With the aid of a more elaborate calculation, the 


difference for y, is lowered to 6 2/60 but yy remains practically unchanged. 


However, as far as we 


know, these results are the best obtained up till now. 


Les calculs relatifs aux susceptibilités magné- 
tiques principales du benzène présentent un double 
intérêt : 

a) La comparaison des résultats théoriques et 
expérimentaux est 1ci possible non seulement pour 
la susceptibilité moyenne, qui était la seule gran- 
deur mesurée dans le cas des molécules envisagées 
jusqu'à maintenant [1], mais aussi pour l’aniso- 
tropie qui a pu être atteinte par des mesures sur 
des monocristaux de benzène [2]. Ceci nous permet 
de mieux apprécier les conséquences des approxi- 
mations nécessitées par la théorie [3]. 

b) L’anisotropie calculée du benzène peut servir 
de référence pour la détermination en valeur 
absolue de l’anisotropie d’autres molécules de 
nature polycyclique, puisque nous disposons de la 
méthode de London qui conduit seulement à des 
valeurs relatives rapportées à l’anisotropie. du 
noyau benzénique [4]. Il convient de noter que 
divers essais ont déjà été effectués pour connaître 
l’anisotropie absolue théorique du benzène, soit 
par la méthode de Pauling qui cherche à adapter la 
formule classique de Langevin [5], soit par des 
modifications de la méthode originale de 
London [6]; ceux-ci dépendent toutefois d’un 
cerbain nombre d’hypothèses simplificatrices, en 
particulier de celle qui attribue l’anisotropie aux 
seuls électrons x en négligeant le rôle des élec- 
trons 6 (les résultats sont généralement trop faibles 
dans ces conditions). 

Nous avons donc entrepris le calcul des trois 
susceptibilités principales de la molécule de ben- 
zène en nous servant de la théorie générale [3] qui 


nous avait déjà fourni de bons résultats dans de 
nombreux cas [1], et nous supposons que la fonc- 
tion propre de la molécule (dans son état fonda- 
mental) est correctement décrite par un produit 
non antisymétrisé d’orbitales monoélectroniques, 
elles-mêmes combinaisons linéaires d’orbitales ato- 


miques (méthode des orbitales moléculaires 
L C A O). 
I. Susceptibilités des électrons s et o. — Les 


calculs relatifs aux susceptibilités des électrons 1s 
du carbone (on suppose que ceux-ci ne participent 
aucunement aux liaisons) ou à celles des liai- 
sons & CG —C et C — H ne présentent pas de parti- 
cularités à signaler par rapport aux déterminations 
déjà effectuées dans des cas semblables [7]; il 
convient seulement de tenir compte de la variation 
des distances internucléaires. 

Les orbitales moléculaires localisées décrivant 
les liaisons © dériveront des orbitales atomiques 
hybrides trigonales du carbone tandis que les dis- 
tances choisies seront les suivantes : : 


CE 1539 À 0 et CPE TOR 


Quant aux fonctions orbitales atomiques, leurs 
formes analytiques sont celles de Slater [8], avec 
les charges nucléaires effectives indiquées par 


.Coulson [9]. Enfin, pour les liaisons C — H dont la 


polarité est mal connue, nous admettons que l’orbi- 
tale moléculaire localisée 


Ÿ = N(ÿc + Vr) 


est acceptable, où dc et Lx sont respectivement les 


_ N°7 


orbitales atomiques normalisées convenables asso- 
ciées au carbone et à l’hydrogène. Les résultats 
obtenus avec ces hypothèses se trouvent réunis 
dans le tableau T où 11 et y, désignent les suscep- 
tibilités correspondant à la direction parallèle et 
aux directions perpendiculaires à l’axe de la 
haison. 


TABLEAU I 


(— xm.10$ u. e. m.) 


Lraisons : L + E 

(A 2 ÉLECTRONS) Xi x Lx 4 x 
+ Électrons 1s de C 0,45 0,15 015 0,15 
CC 2 OVNS SN 0 78 0 3,05 12.78 
CH DA ES OA 01 630€ 
IT. Susceptibilités des électrons x. — Si l’on 


note par di (i — À, B, C, D, E, F) l’orbitale ato- 
mique 2p, (cf. fig. 1) attachée au carbone i, les 


ÿ 


trois paires d'électrons x du benzène peuvent être 
décrites par l’ensemble suivant [10] : 

— une orbitale moléculaire totalement symé- 
trique, non dégénérée 


da = Nilba + ŸB + Vo + Ÿp + ŸE + Vr) (1) 


— deux orbitales dégénérées, par exemple 


Ve N(Ÿa in ŸB — VD VE) (2) 
da = Nalba — ŸB — 2 %o — ŸD + Ve + 20r). (3) 


A. SUSCEPTIBILITÉS POUR L'ORBITALE SYMÉ- 
TRIQUE. — Les termes de Langevin 


] 


L e? 


lu — 5 fu + war + (4) 
2 
fun) 1 


ont tout d’abord été évalués en tenant compte de 
l'existence d’un plan nodal pour W,, coïncidant 
avec le plan des divers noyaux, car il convient de 
conduire séparément les calculs dans chacun des 
sous-espaces limités par des surfaces nodales en 


L 
(a, v,w = 29,5; TRES 
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choisissant pour chacun de ceux-ci une origine 
située au centre de gravité électronique partiel [3]. 
Les termes correctifs paramagnétiques définis par 
Liu = Juu — ZE, sont ensuite déduits de la formule 


L L\2 

NT (ko — Xw) 

sr rare 7 L (5) 
Xo + Xw 


déjà utilisée pour toutes les liaisons 6. On remar- 

quera que la symétrie Dyn conservée par , 
A TNT L LES 

entraine y y d'où). —0) 


B. SUSGEPTIBILITÉS POUR LES ORBITALES DÉGÉ- 
NÉRÉES. — [application du champ magnétique H 
a pour effet de lever la dégénérescence existant 
entre les orbitales moléculaires 4, et 3. À partir 
de celles-ci, il convient donc de former des combi- 
naisons linéaires adaptées à la perturbation, c’est-à- 
dire telles que la matrice associée au potentiel 
perturbateur (potentiel comportant un terme en H 
et un terme en H @ H) soit diagonale. Deux cas 
doivent maintenant être distingués, suivant que H 
est appliqué normalement ou parallèlement au plan 
de la molécule. 

1) Champ H normal au plan du cycle benzénique. 
— Il est assez facile de vérifier, par de simples 
considérations de symétrie, que 


_<hMih? #0() 


. k d d Jtsuf + 
OÙ, = Di Ê > y à est l’opérateur associé 


à la composante suivante Oz du moment orbital. 
On a, d’autre part 


< Vol M]Ÿe PIE Yale] Ÿs EN) 


et la dégénérescence se trouve levée par la pertur- 
bation du premier degré en H. Les fonctions adap- 
tées à la perturbation seront nécessairement 


da = 72 (de £ id). 


Ces fonctions présentent l’inconvénient de deve- 
nir complexes alors que la théorie relative au 
calcul des termes correctifs y+ a été établie seu- 
lement pour des fonctions réelles pures. Il est utile 
de noter que ces deux nouvelles orbitales molé- 
culaires correspondent à la même densité élec- 
tronique puisque 


DE Ye = VE = Ye dl + VE (6) 


De plus, cette densité respecte la symétrie Don 
comme on peut le voir directement à partir des 
relations (2) et (3) après calcul des coefficients de 


normalisation V, et N; (V3N, — N;), ou encore 


G) Notations : < Y,[A|W > = fui A ga, etc. 
Lan 
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en vérifiant que 4 et d_ admettent les expressions 
équivalentes 


dr = Naulba + e—ir8 bp + e—°i7/8 Vo 


L e—ir Yn + e— 57/8 Qn + e—Sir8 y] 
DENN = E PPRee mEnbe 
L er Qp + ets Qu + esir/8 Vr] (7) 


AvOC Ne = IN, 

D'autre part, chacune des orbitales D, et 4 
considérée séparément donne lieu à un moment 
magnétique permanent mais les effets de ceux-ci 
se compensent exactement, par suite de légalité 


< G+lMlbr > =— <I> #0 (8) 


et étant donné que les orbitales se trouvent comple- 
tement occupées : le moment permanent résultant 
continue par conséquent de rester nul. 

Après ces quelques remarques, il devient possible 
de poursuivre la méthode de calcul du coefficient 
correctif y* Pour ce cas particulier, si nous prenons 
comme précédemment une fonction perturbée en 
présence du champ %, identique à d, multipliée 
par (1 + a.X,.xy), à cette différence près que le 
paramètre a ne sera plus nécessairement imaginaire 
pur. On trouve dans ces conditions (pour deux 
électrons) | 


2 
22 ni < dalx? + y?|h4 > 
2 
se [a* < delay M,] be > + à < dilazyM,| + >] 
h?|a|? 2 
NT = _ « d+|r? Me y?| + 20 (9) 


Or, étant donnés les caractères de symétrie des 
fonctions 4 et 4, le terme entre crochets de 
lexpression (9) (du premier degré en a) est nul, si 
bien que l’extremum de 7x sera obtenu pour a — 0. 

Dans le cadre de notre approximation actuelle, 
nous n'avons pas à tenir compte ici d’un terme 
correctif (exactement comme dans le cas de l’orbi- 
tale symétrique) et le calcul se réduit finalement à 
clurde y =". 

2) Champ H parallèle au plan du cycle benzénique. 
— On peut tout d’abord vérifier sans grandes diffi- 
cultés la nullité des éléments matriciels suivants : 
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ce qui signifie que la dégénérescence ne se trouve 
pas levée par les termes du premier degré en H 
appartenant au potentiel perturbateur. On. a, 


d’autre part : 
<< ble? + 2h3 > = < Yely? + ds > =0 
tandis que les éléments 
< Yale? + 212 7, 

< plz? + 22|ds >, 


seront tous différents de zéro. 

Nous constatons de cette manière que la dégéné- 
rescence est effectivement levée par les termes 
perturbateurs en H @ H;et ce sont précisément 
les fonctions %, et d, primitivement choisies qui 
conviennent pour la direction du champ H actuel- 
lement étudiée. 

En dehors du plan xOy, les fonctions 4, etd, 
possèdent chacune un autre plan nodal perpen- 
diculaire au eyele : ce plan passe par les sommets C 
et F pour », ;ilest médiateur des liaisons AB et ED 
pour %, (cf. fig. 1). Il n’a malheureusement pas 
été possible de tenir compte de ces plans nodaux 
d’une manière totalement rigoureuse et nous 
n’avons pu déterminer la position du centre de 
gravité électronique de chaque espace partiel que 
moyennant une légère erreur par défaut. Cette 
difficulté numérique introduit une petite erreur 
par excès au cours de la détermination des suscep- 
tibilités principales. 

Par suite de la symétrie sénaire présentée par la 
densité électronique globale associée aux électrons 
occupant les orbitales 4, et 4, la susceptibilité 
correspondant à cet ensemble de 4 électrons doit 
être la même quelle que soit la direction considérée 
pour H, dès que celle-ci appartient bien au plan 
de la molécule. Cette remarque impose l’égalité 


Xaa( Ve) + XV) = Xyy( Ve) + Xuyl Vs) 


qui serait rigoureuse si %, et V, étaient les fonc- 
tions propres exactes adaptées au problème et si les 
formules de calcul des susceptibilités n’étaient pas 
approximatives. Comme ce n’est pas le cas, nous 
avons entrepris séparément le calcul avec les deux 
fonctions et nous avons utilisé ultérieurement 
comme évaluation de Yzx — Yyv là plus petite des 
deux valeurs absolues obtenues et correspondant à 


< ph + M >, 
< plu + 343 >, 


(10) 


< 4|M. AT 
VelMalgs 7, < elMlhs 7, l’un des membres de l’égalité (10). Cette manière 
< Ye Mol be >, < YMilhs >, de procéder est conforme aux principes généraux 
TABLEAU II 
(— x.106 u. e. m.) 
ORBITALE 
(complète à 2 électrons) Lee XL x xE Le a Xxx Xuv 
Yi(totalement symétrique) 17,36 0,78 006.68 NE 2,87 2,87 
1,34 4,94 2,56 #0 
dégénérées [a 19,07 1,04 DR 
8 4,41 D TRE 020 — 9,10 2,79 3,35 
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de la méthode de variation ; elle doit également 
permettre de corriger dans une certaine mesure 
l'effet de l'erreur par excès précédemment signalée. 
Les divers résultats obtenus se trouvent indiqués 
dans le tableau IT. 

On notera que c’est le total 


Xan(% 2) de Xuu(Ÿ3) = 5,09.107$ 


(cf. dernière colonne du tableau IT) qu’il convient 
de retenir pour la susceptibilité des 4 électrons 
occupant les orbitales dégénérées lorsque H est 
parallèle au plan de la molécule. Il est par ailleurs 
intéressant de constater que cette valeur est très 
prophende="5,74 1078 — 2;%,,(0:) — 2y(0). 
Avec la technique d’approximation actuellement 
utilisée, les contributions de %,, d,, et 4, à la 
susceptibilité y, (4 — x ou y) nous apparaissent 
comme très voisines. 


IIT. Susceptibilités principales, susceptibilité 
moyenne et anisotropie du benzène. — Les résul- 
tats groupés dans les tableaux I et IT nous per- 
mettent maintenant d’obtenir les susceptibilités 
principales de la molécule de benzène, en tenant 
compte de la nature tensorielle de ces grandeurs. 
Nous avons pour la direction normale au plan du 


| cycle : 


Xe = 6x(1 C) + 6x1(CC) + 6yL(CH) + Xl di) 
+ 2X#(V+ ou bd) 
et 00-110%e.. 


(11) 


Ensuite, pour toute direction parallèle au plan 
du cycle, nous trouvons : 


Xax — Xuy — 6x(1s C) + 2yn(CC) + 4£xn(CC) (12) 


Cons + 47, (CC) sin? > + 9yn(CH) L-4yn(CH) 


cos? © + 4y1(CH) sin?> + Xyl di) + Xwv( Va) + Xl Va) 


soit numériquement 
étre DT TETE AID 


La susceptibilité moyenne molaire devient dans 
ces conditions : 


1 
AM 3 (Xe ce Dr) TT 65,15.10—$ (13) 


_etl’anisotropie : 


À = fre nl = AO (14) 


Les valeurs expérimentales de comparaison dont 
nous disposons [2] sont yf — — 93,1.107$ et 
4 — — 35,70.1075$, si bien que les valeurs qui 
viennent d’être calculées apparaissent comme enta- 
chées d’une erreur de 7,4 % pour x et de 33,7 % 
POUT Yx — Xw. Ce dernier écart est très nettement 
supérieur à tous Ceux que nous avions constatés 
jusqu’à maintenant et ceci nous a incités à recher- 
cher des améliorations de notre méthode d’approxi- 
mation. - 
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IV. Recherche d’une meilleure approximation 
pour les susceptibilités liées à l’orbitale x tota- 
lement symétrique. — Étant donné que les résul- 
tats numériques ont toujours été satisfaisants pour 
les composés ne comportant que des orbitales de 
type o, 1l paraît vraisemblable que la plus grande 
cause d'erreur provient de l’estimation des suscep- 
tibilités associées aux électrons x. Un examen, 
même superficiel, de la forme analytique de %, 
d’une part, de 4, et 4, (ou de d+ et d_) d'autre 
part, plaide en faveur de cette hypothèse car la 
répartition des électrons ne présente plus une 
symétrie de révolution autour d’un axe, mais seule- 
ment une symétrie sénaire par rapport à l’axe 
passant par le centre de l’hexagone benzènique et 
normal à son plan. 

Pour la recherche d’une meilleure technique 
d’approximation, nous nous sommes bornés à 
l’étude de l’orbitale b.,, celle-ci présentant l’avan- 
tage d’être réelle et de ne pas introduire de compli- 
cations dues à la dégénérescence. Les compor- 
tements respectifs de 4, et des orbitales dégénérées 
au cours de la précédente évaluation (cf. tableau II 
et fin du paragraphe IT) nous laissent en effet 
présager que leurs contributions resteront du même 
ordre de grandeur, ce qui devrait conduire à des 
corrections relatives analogues. Pour ce nouveau 
calcul concernant d.,, nous avons utilisé une forme 
modifiée de la théorie générale, forme mise au point 
à propos d’une étude critique des approximations 
utilisées jusqu'à maintenant [11] et dont nous 
rappelons très brièvement les grandes lignes. 

En introduisant la fonction (connue en prin- 
cipe) F,(o,) ne dépendant que de la variable angu- 
laire o, (angle de rotation autour de l’axe Ou) et 
définie par l’intégrale 


CO +00 3 ,2 
Fy(@u) £ he 5 4? du du 


(les cocrdonnées utilisées 4, p4, ®, sont de type 
cylindrique avec Ou comme axe polaire) 

eten admettant que la fonction G,, définie par l’orbi- 
tale perturbée %,[1 + iJ,.G,] en présence du 
champ %,, est de la forme 


(15) 


Guy — pu Du Eu) (16) 


il apparait que la nouvelle fonction inconnue ®,(o,) 

doit satisfaire à l’équation différentielle suivante 

obtenue par application du principe de variation : 
d? ®, . d(Log Fy) d® d (Log Fy) 
doi dPy dou dou 


| 
Ï 


. (17) 


& Dy — 


C'est cette équation (17) que nous allons main- 
tenant utiliser de deux manières différentes suivant 
la direction supposée du champ H, après avoir 
remarqué que l’approximation employée jusqu’à 
présent consiste à admettre que 


DalPu) = à Sin 2Pu (138) 


est acceptable comme solution approchée de (17). 
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À. SUSCEPTIBILITÉ SUIVANT UNE DIRECTION NOR- 
MALE AU PLAN DE LA MOLÉCULE. — Nous avons 
vu au paragraphe II-A que le terme correctif para- 
magnétique yÀ, est égal à zéro lorsque l’on s’en 
tient à l’approximation (18). Par suite du rempla- 
cement des axes d’ordre 2 (présents dans les 
liaisons o homonucléaires) par un axe d’ordre 6, la 
transformation la plus naturelle de l'hypothèse (18) 
consiste à poser 


®D,(o.) — «à sin 6. (19) 


Cette dernière relation étant admise, l’équa- 
tion (17) pourra réciproquement servir à déter- 
miner F,(v.), si bien que l’on arrivera à se faire 
une idée de la validité de (19) en comparant (au 
moins pour certaines valeurs de l’angle ®.) les deux 
fonctions F, obtenues, l’une à partir de la défi- 
nition (15), l’autre par l’intermédiaire de (17). De 
cette dernière équation, présentée sous la forme 


d (Log F;) __  h0xsin 6% 
dœ, 7 1 — 6x cos 69, 


10 d(1 — 6x cos 6;) 


9 ° 1 — 6 cos 6% 


on tire immédiatement 


K 
(1 — 6 cos 6w,)10/9° 


F( 92) = (20) 


Cette forme analytique est bien susceptible de 
rendre compte des variations assez faibles de F, 
entre ses extrema successifs qui doivent être écartés 
de 300 en respectant la symétrie du noyau benzè- 
nique. Quant aux deux constantes Æ et æ&, nous 
pouvons les obtenir à l’aide de (15) en calculant la 
valeur de F, pour un maximum et un minimum 
successifs, soit pour @, = 0 et w. = x/6 : cepen- 
dant, 1l nous sera suffisant de déterminer unique- 
ment « par l’intermédiaire de la relation 


É s he A0) 

160)  F,{r]6) 
a) Calcul de F,(0). — Le point courant M se 

trouve constamment dans le plan yOz passant 


par les deux atomes de carbone À et D et perpendi- 
culaire au noyau benzènique (plan de la figure 2). 


(21) 


Ce sont ici presque uniquement les orbitales 
atomiques da et Yn qui donnent la densité élec- 
tronique et nous négligerons les contributions 
faibles de Vs, Vo, br et dr. Avec les notations de 
la figure 2, nous obtenons d’après (15) 


F,0) = M [°° [UT eilta + vo}? des àz (22) 


CO Me 
#2 [7 dre f 'rtexp(— 2e V2 + (oz — R)°] dz. 


L'intégration par rapport à z est analytiquement 
possible et conduit à un premier résultat s’expri- 
mant à l’aide des fonctions de Hankel d’ordres 0 
et 1 et d’argument 2ico, [12]. Quant à la deuxième 
intégration par rapport à p,, elle a été effectuée 
numériquement. 

b) Calcul de F(x/6). — Le point courant M est 
cette fois dans le plan médiateur de AB, et nous 
avons pris comme forme approchée de la fonction 
propre 


NilŸa =F ŸB) Re 2N; Ya. 


En utilisant les notations indiquées sur la fi- 
gure 3, on a 


ha AVS) IE 
TA = TB = [+ #4 es) 


rG08: 


et l’intégrale de définition (15) devient 


+ 
Erl6) = sx? [ EU 


+0co : [ — 
{ 3? exp L— 2c V/:? + p2— PÈRV3 + R?] az. (23) 


N°07 


NY 


. Le procédé de détermination de (23) est iden- 
tique à celui indiqué ci-dessus pour (22). 
©) Valeur numérique de «. — Les calculs qui 
viennent d’être schématisés donnent pour le rapport 
de F,(0) à F,(x/6) la valeur 2,045 et la résolution 
de l’équation (21) nous permet d’évaluer x — 0,052. 
d) Calcul de la susceptibilité y. — Une fois le 
coefficient « connu, nous pouvons évaluer la suscep- 
übilité x de la manière suivante en fonction du 
seul terme de Langevin y£. On a en effet [11] 


e? ET d®, 
af e(—a)n 6 


soit encore, en faisant apparaître y£ : 


ÉLIETE 


DE C7 © 2 - (24°) 


Lorsque nous utilisons  l’approximation 
D, — « sin 6,, nous arrivons de cette façon à 


1 ut) d œ. 
L JO (1 — 6% cos 69,)1® 
Ka =. Jia 5 


(25) 


T/6 do, 
ÿ (1 — 60 cos 69,)10/9 


La correction définitive étant faible, nous pou- 
vons nous permettre d'évaluer graphiquement les 
deux intégrales présentes dans (25), ce qui nous 
donne 


CT 6,355 


CN ES race 


Tandis que l’étude effectuée au paragraphe IT-A 
nous donnait automatiquement un coefficient cor- 
rectif paramagnétique y£ égal à 0, nous avons 
cette fois-ci 


L € L 
à —\ #29 Xax — 0,118 y:7, 


si bien que la susceptibilité calculée y£ relative à 
l’orbitale totalement symétrique 4, prend main- 
“tenant la valeur — 15,40.1076 au lieu de 
— 17,36.10—$ (cf. tableau IT). 


B. SUSCEPTIBILITÉ SUIVANT UNE DIRECTION 
PARALLÈLE AU PLAN DE LA MOLÉCULE. — (C’est 
dans cette direction que l’écart vis-à-vis des résul- 
tats expérimentaux est le plus considérabla et 
nous avons tenté la résolution complète de 
l'équation (17) après évaluation numérique de la 
fonction (15). Cette partie de notre travail a été 
possible grâce à l’Institut Européen de Calcul 
Scientifique qui a bien voulu mettre à notre dispo- 
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sition un certain temps de calcul sur l’Ordi- 
nateur IBM 704. Nous tenons à remercier ici cet 
Institut pour son extrême obligeance. 

Lorsque le champ H est appliqué parallèlement 
à la direction Oy (cf. fig. 1), nous allons utiliser 
l'équation (17) sous la forme : 


e Le, (i +) + 4Fi(eu). Di(pu) = 0. (26) 


y 


Par intégration des différents termes de (26) 
entre les valeurs 0 et ,, nous aboutissons à l’équa- 
tion intégro-différentielle équivalente permettant 
de déterminer ®, : 


d®,\1? p 
[Fulen) (1 ml He ‘ y D Ar — 0 


y 


(26°) 


Cette dernière équation présente l’avantage de 
permettre le calcul numérique à partir des seules 
valeurs de F,(v,) en évitant la détermination préa- 
lable de dÆ,/dw,, et sa résolution a ensuite été 
conduite à l’aide des étapes ci-après : 

a) Calcul de la fonction Fy(o,). — Il convient 
d'évaluer tout d’abord l'intégrale (15) pour des 
valeurs de +, suffisamment rapprochées et compa- 
tibles avec les pas d’intégration que nécessitera 
ultérieurement la résolution de (26). Or, il apparaît 
d’assez sérieuses difficultés analytiques qui pro- 
viennent essentiellement des termes rectangles du 
type Va YB, Ya Vo, etc. présents dans le dévelop- 
pement de df. Cependant, l’évaluation finale 
de yyy au moyen de la formule (24) (avec substi- 
tution des indices y aux indices z) nous montre 
que F, pourrait sans inconvénient être remplacée 
par n'importe queile autre fonction qui lui soit 
proportionnelle. La forme de (15) nous a donc 
incités à négliger tous les termes de recouvrement 
dans l’expression de la densité électronique car la 
fonction qui se déduira de cette simplification 
variera comme la fonction F, exacte, en lui restant 
d’ailleurs constamment inférieure. Par suite, pour 
ce calcul, nous utiliserons 


DA NA VA + VB + dé + 4 + VE + VE). 


Pour tenir compte du fait que le plan xOy est 
un plan nodal, il convient de mesurer la coor- 
donnée z à partir du centre de gravité électro- 
nique G du demi-espace supérieur D* par exemple 
(z > 0), défini par la relation donnant sa cote soit 


za 


Nous devons maintenant calculer l'intégrale (15) 
dans un plan MGy’ parallèle à Oy et coupant xOy 
sous l’angle o, (cf. fig. 4). Cette intégrale va 
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Trace du plan 


Free 


pouvoir être décomposée en une somme de 6 inté- 
grales monocentriques, relatives respectivement 
aux sommets À et D,B et C, E et F. Nous aurons, 
avec les notations de la figure 4 : 


Res la CR CRT Or DEC) 


en posant par exemple (le coefficient V2 pourra 


être omis sans inconvénient d'après les précé- 
dentes remarques) 


+co +0co 3 
= lie eu vA dey dy 
+co f+co À 
=) jh AZ + pu sin 2 
0 —co 


exp [— 2e V(Z + po, sin 9)? + pè cos? , + y?] de, dy 


Fy D (pu) | = Fyl Py) 


(28) 
lorsque l’angle &, est compris entre 0 et x/2 
—%/sinp +co 
File) = F6) = [7 [7e dA de, dy 
pour — rf2 < p, < 0. (29) 


Les intégrales doubles précédentes sont réduc- 
tibles à une seule quadrature. En posant 


u—=Z +eysin ou —ulte y et Z' = Z]tg o, 
on trouve pour (28) 


T HOME 75 
Her = u?{u — Z) 


HT [2ie Vur + (uw — 7j] du (28’) 


Fi (y) = Var + 27}: 


tandis que le même type d'expression est valable 
pour (29), les bornes d'intégration devenant Z 
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et 2Z au lieu de Z et + co. Dans toutes ces for- 
mules, H® représente la fonction de Hankel 
d'ordre 1. 

Les fonctions FE — F£ se trouvent encore 
données par les mêmes quadratures à condition de 
prendre Z’ = R VD TO RCRULE oo 
la longueur des liaisons GC — C du noyau). Quant 
aux fonctions FF et FF, leur valeur commune 
est tout à fait négligeable devant les autres dans 
l'intervalle — x/2 < o, < + x/2: 

Il convient de noter que les quadratures du 
type (28) ne sont plus utilisables pour ®, = 0. 
L’indétermination peut heureusement être levée 
grâce aux relations valables dans ce cas particulier. 


FAO) = — xZ? es WVZS Eu Hoi VZ? TE ui] du 


(30) 


(u— RV3/2) 


"ù uÿ V/2? eo 


HD oie V/22 + 


CV) du. (30) 


L'évaluation numérique des intégrales (28), (30) 
et (30/) n’a pu être effectuée qu’à l’aide de l’Ordi- 
nateur IBM 704. Nous avons appliqué la méthode 
de Simpson [13] et déterminé F,(w,) de 1/2 degré 
en 1/2 degré ; la figure 5 correspond à la courbe 
ainsi obtenue. Quelques difficultés pratiques sont 
relatives au choix du pas d'intégration et c’est 
malheureusement au voisinage de son maximum 
que la courbe est la moins précise ; quelques irrégu- 
larités apparaissent également dans cette région 
par suite du passage nécessaire de la formule (28') 
aux formules (30) et (30). 


FYP,) 


Lu IT 


Nes 6 


YA 2 LM 


b) Obtention de la fonction inconnue ®D,(o). — 
La résolution de l'équation intégro-différen- 


NO) 


telle (26) a également été entreprise numéri- 
quement. La principale difficulté concerne les con- 
ditions initiales car nous devons imposer à d, 
d’être périodique avec T — 2x. F, étant une 
fonction paire, les considérations de symétrie nous 
indiquent bien que ©, doit être impaire d’où 
D,(0) — O0 mais nous ignorons complètement la 
valeur de (0). 

Par suite, il convient de calculer deux solutions 
particulières ®, et ®, répondant à 


®,(0) Sy D,(0) = 0, 


la périodicité devant ensuite être imposée à la 
combinaison linéaire ®, — a.®, + b.®,, les para- 
mètres a et b se trouvant liés par «a + b — 1. En 
pratique, les solutions particulières ont été cal- 
culées par la méthode de Runge-Kutta [14], à 
l’aide du calculateur électronique Gamma de la 
Compagnie Bull (*) et le problème a été poursuivi 
à l’aide de l’Ordinateur IBM 704 de la Société IBM 
France (*) selon la méthode d’Adams. Les diff- 
cultés relatives aux pas d'intégration ne nous ont 
toutefois pas permis d’arriver à une solution entiè- 
rement satisfaisante car nous avons constaté que 
la dérivée seconde ®,(w,) n’était pas exactement 
nulle pour +, — 0 ni surtout pour ®, = tx, en 
même temps qu'une certaine irrégularité (aplatis- 
sement de la courbe) apparaît pour ®'(x/2) 
(voir fig. 6 et 7). 

c) Calcul de la susceptibilité yyy. — Connaissant 
le terme de Langevin y£, (cf. tableau IT), nous 
maintenant obtenir la 


pouvons susceptibilité 
suivant l’axe Oy par 
me M 
7e Fuleu)lt — D] dey 
L Jo 
Xuu.— Xuv: = 
. Fyleu) dy 
0 
09/46 10m 00/3222 5-05 105 


Ce résultat peut paraître en contradiction avec le 
théorème fondamental de la méthode de variation 
puisque nous devrions obtenir une valeur absolue 
inférieure à 2,87.107$ (cf. tableau IT) en intégrant 
directement l’équation (26’). En réalité, cette diver- 
gence n’est qu’apparente et due seulement au fait 
que la fonction ®, utilisée est seulement très 
proche de la solution rigoureuse, ainsi que nous 
avons pu le constater en effectuant de très légères 
modifications empiriques de ®, et en employant 
la formule complète 


bout sf ol — où + 403] de, (31) 


(*) Qu'il nous soit permis de remercier ici les ingénieurs 
de la Société IBM France et de la Société Bull pour l’atta- 
chement qu'ils ont apporté à notre problème et les nom- 
breux conseils qu’ils n’ont cessé de nous prodiguer. 
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où il n’est plus nécessaire de supposer que .®, est 
solution exacte de (26). Nous avons trouvé de 
cette façon des résultats pratiquement identiques 
à — 2,87.107$ mais cependant sans aucune amélio- 
ration significative. A titre indicatif, pour les fonc- 
tions ®, et ®, qui nous sont apparues comme les 
meilleures, un certain nombre de valeurs numé- 
riques sont portées sur le tableau III. 


TABLEAU III 


Py Fy D, D 
0° 0,0057 0 — 0,89 
11 0,0061 — 0,151 — 0,83 
200 0,0072 — 0,285 — 0,70 
300 0,0091 — 0,390 — 0,51 
400 0,0140 — 0,455 — 0,25 
D 09 0,0261 — 0,469 + 0,07 
60° 0,0746 — 0,421 + 0,49 
2100 0,3309 — 0,315 + 0,71 
80° 0,8906 — 0,182 + 0,80 
90° 1,4345 — 0,039 + 0,83 
1000 0,8908 + 0,103 + 0,80 
111400 0,4460 + 0,231 + 0,66 
1200 0,2076 + 0,329 + 0,46 
1800 0,0998 + 0,387 . + 0,22 
1400 0,0540 + 0,399 — 0,09 
1500 0,0344 + 0,399 — 0,38 
1609 0,0258 +. 0,271 — 0,63 
17100 0,0219 +- 0,144 — 0,79 
1800 0,0202 0 — 0,84: 


F16G.16: 


V. Conclusions. — On voit que si les compli- 
cations de méthode ont très sensiblement amélioré 
les résultats numériques pour la direction perpen- 
diculaire au plan du cycle benzènique, il n°y a par 
contre pas de modification pratique pour les direc- 
tions parallèles à ce plan. Si nous admettons que la 
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correction pour les orbitales x dégénérées est du 
même ordre de grandeur que celle valable pour 
l’orbitale symétrique, ce qui semble logique puisque 
les densités électroniques partielles ont les mêmes 
symétries, nous obtenons en définitive 


Xua = — 93,64 .10—6 (erreur de 0,6 %) 
La = Xyy = — 47,74.10—$ (erreur de 33,7 %) 
x = — 63,04.10—5 et A — 45,90.1076 


Le résultat pour l’anisotropie A peut sans doute 
être considéré comme très légèrement amélioré par 
rapport à celui obtenu par Itoh, Ohno et Yoshizumi 
à l’aide de la méthode de London [6] pour les 
électrons x, soit À, — 38,2.10756 car en ajoutant 
à ce À, 5,28.1076, qui correspond suivant nos 
calculs à la contribution des électrons © à l’aniso- 
tropie, on trouve seulement À = 43,5.107$ au 
lieu de 57,4.1076 expérimental. 


Manuscrit reçu le 8 avril 1960. 
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PHÉNOMÈNES DE RELAXATION DIPOLAIRE EN FONCTION DE LA TEMPÉRATURE 


Par S. LE MONTAGNER, 


Faculté des Sciences de Rennes. 


Résumé. — Les équations de Debye de relaxation dipolaire sont examinées en prenant la 
température comme variable indépendante. Les prédictions théoriques donnent des formes de 
courbes e(T) sensiblement différentes de celles e(w), notamment dans les régions de température 
où la dérivée d:,/0T est négative, c’est-à-dire au voisinage d’un changement de phase. Discussion 
et vérification expérimentale sur deux aluns d’ammonium. 


Abstract. — The Debye dipolar relaxation equations are examined taking temperature as the 


variable. 


Introduction. — Le phénomène de ‘relaxation 
dipolaire dans les diélectriques solides ou liquides 
soumis à un champ électrique alternatif est bien 
connu et il est en général convenablement repré- 
senté par les équations de Debye donnant les 
parties réelles et imaginaires €, et e, de la cons- 
tante diélectrique complexe e*. 


Est Eos 
L£1 — £oco Lio 73 d2 mr (1) 
HE OT : : 


en fonction de la constante diélectrique statique e, 


(à fréquence nulle) ; du carré de l'indice de réfrac- 
tion nr? — <, de la fréquence angulaire © et 
d’une quantité + appelée temps de relaxation que, 
pour la commodité de l’exposé, nous supposerons 
ici unique et obéissant à la loi : 


r = Toexp (UJXT) (3) 


ro €t U énergie d’activation étant des constantes 
dépendant du corps considéré [1], [2], [5]. 


Br 


e, et e, sont donc fonction des deux variables w 
et retilest, pour le calcul, beaucoup plus commode 
de prendre © pour variable, T pour paramètre, 
1.e de représenter à température fixe les variations 
de €, et e, en fonction de la fréquence, d’où les 
courbes classiques (fig. 1). 


1 >. Theoretical curves e(T) are predicted, somewhat different from the classical one e(w) 
chiefly in the regions where De,/07 < 0, i.e near a phase transition. 
were used to make an experimental verification. 


: } 


Two ammonium alums 


Expérimentalement, toutefois le problème est 
tout autre, et il est très difficile de pouvoir faire 
des mesures sur un très grand intervalle de fré- 
quences, alors qu’il est commode de faire varier la 
température à fréquence fixe. Cette solution est 
adoptée par de nombreux auteurs, mais il ne 
semble pas que, parallèlement, une investigation 
ait été menée pour chercher les prédictions corres- 
pondantes des équations de Debye. C’est ce que 
nous allons nous efforcer de faire, en supposant 
explicitement que dans l’intervalle de température 
exploré (sauf peut-être aux bornes de celui-c1), le 
corps étudié ne présente pas de changement de 
phase. 


I. Étude de €, à fréquence fixe en fonction de 
la température. — 1.1. POSSIBILITÉ D'UN MAMXIMU 
DE LA COURBE €,( 7°). — La forme de la courbe e.(T') 
à fréquence constante doit être a priori, moins 
simple que celle de e;(w) à température fixe, 
puisque les quantités «,, & et + sont fonction de la 
température seule. 

Nous savons toutefois que & varie généralement 
peu avec la température, et nous le considérerons 
ici comme indépendant de 7. Posons & — ss = 4 
et æ' — oxfoT, quel que soit x. 

La dérivée de,/dT peut se mettre sous la forme : 


[8] 
= [ét +ots)+ et] (0) 


si e, est positif, #i l’est également et la variation 

de e.(T) est monotone, constamment croissante, 

rappelant à une inversion près, la courbe e;(w) 
précédemment rencontrée. 

Si par contre &, est négatif e; peut s’annuler à la 

température T', racine de l’équation implicite en T': 

2 2 
es (1 + w?7r?) + — = (i}. (5) 


A la température T, il y aura un maximum de la 
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constante diélectrique réelle &,(T), ce qui constitue 
une différence essentielle avec les courbes e,(w). 


4.2. VARIATION AVEC LA FRÉQUENCE DE LA POSI- 
TION T, DU MAXIMUM DE €,(7”). — L’équation (5) 
transformée par (3) s'écrit : 


w? 2 exp (2UJRT) [e’ + (2aUJKT?)] + & = 0. 


Dans un petit intervalle de température, on peut 
toujours supposer une variation linéaire pour e;(T) 
d’où e, — Cte. Développant en série l’exponen- 
tielle, on voit que le coefficient de w? est une 
fonction décroissante de la température, d’où néces- 
sité que le sens de variation de T', soit le même 
que celui de ©. 

Donc la température T, du maximum de e, est 
une fonction croissante de la fréquence. 


1.3. CALCUL DU MAXIMUM DE €; — Notons 
le Eimax. 11 se calcule en éliminant wr entre les 
équations (5) et (1), soit 

é Em = © + (AT? J2U) &. (6) 


Comme l’on suppose ici nécessairement & <0 
(pour qu’il y ait un maximum) e,, est inférieur à &, 
à la même température, et la courbe e,(T) relative 
à une fréquence donnée se trouve en dessous de la 
courbe e(T), au moins au voisinage de son 
maximum. 

On peut montrer du reste qu'il en est ainsi pour 
le reste de la courbe e,(7), en effet pour un corps 
polaire on a toujours : 


ÉCOLE 
£oo &° T°? < eg W? T°? 
Eco + @?T?) — eco + € < e5(1 + w°? +?) 
D'où par division par (1 + w? +?) 
CC 


De plus, la différence es, — €,» est comme 7, 
une fonction croissante de la fréquence, et pour 


dE, 


d’où 


Hire 2) 


une fréquence donnée (le long d’une courbe expéri- 
mentale) une fonction décroissante de la tempé- 
rature, tendant vers 0 si T augmente. 

La courbe e,(T) est donc l’enveloppe de toutes 
les courbes €; à partir d’une température suff- 
samment grande devant T.. 
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L'aspect schématique de cette famille de courbes 
est donc celui de la figure 2. 


1.4. Discussion. — L'hypothèse €; < 0 dans 
un domaine de température suppose [5] une aug- 
mentation d’entropie pour le corps étudié lorsque 
la température décroît, c’est-à-dire implique la 
nécessité d’un changement de phase à une tempé- 
rature inférieure, soit T4. 

Selon les valeurs numériques des quantités figu- 
rant dans l’équation (5) la racine T, peut être 
supérieure ou inférieure à 7. Dans le 127 cas, les 
résultats prévus ci-dessus sont valables et on devra 
observer effectivement un maximum de €, à la 
température T,. Si au contraire on a T; <T,, la 
solution est virtuelle, puisque à 7, il y a variation 
de structure, entraînant des discontinuités des 
paramètres € et de leurs dérivées. Le caleul n’est 
plus valable et la cônstante diélectrique €, ne 
passe plus par un maximum. Du reste pour les 
températures inférieures à T,,on a de nouveau en 
général &; > 0et l’on retrouve le cas banal. 

Les courbes expérimentales présentent alors un 
point anguleux (ou une discontinuité finie) à la 
température 7, du changement de phase analogue 
à celui de e,(T) (suivant que l’on se place dans le 
cas d’un changement de phase du second ordre ou 
du premier ordre). 

D'ailleurs T, étant une fonction croissante de la 
fréquence, 1l existera une fréquence critique 
tellé que pour © = ©," on ait To) = "76e 
au-dessus de laquelle les courbes e,(7) présentent 
à nouveau un maximum. 

En résumé, pour un corps présentant un chan- 
gement de phase à la température T, : la constante 
diélectrique statique comme la plupart des para- 
mètres d’état présente à cette température un point 
anguleux ou une discontinuité finie ; la partie 
réelle de la constante diélectrique complexe à une 
fréquence © peut, soit présenter une variation 
entièrement analogue, soit présenter en outre un 
maximum à une température variable avec la fré- 
quence et en tout cas supérieure à T,. 

Expérimentalement, il est important de ne pas 
confondre ce maximum de €; avec un point angu- 
leux ce qui conduirait à une valeur erronée pour la 
température du changement de phase. 


IT. Étude de :, à fréquence fixe, en fonction 
de la température. —— 2.1. TEMPÉRATURE DU MAXI- 
MUM DE £,(T). — Il sera défini par la condition 


dE: JO T— 0 
soit 
& WT 7e at’ 2a a T7’ 
Se A DE TE w?r? {1 2 ASE 0 


| 
| 
| 


% 2 _9 
Es OT EN ire 
nes OT 
1 ET: 


(7) 
fournissant implicitement par (3) la température T, 
du maximum de e,. 

On note immédiatement que n’est plus vérifiée 
la condition classique de relaxation ©+ — 1, pour 
laquelle à la température 7 on observe en fonction 
de la fréquence, un maximum de e, à la fréquence ©. 

Il est donc indispensable pour pouvoir comparer 
des mesures faites en fonction de la fréquence, 
avec des mesures faites en fonction de la tempé- 
rature, de repérer cette température 7. 

Soit une mesure faite en fonction de la fréquence, 
à la température T;, pour laquelle e,(w) a son 
maximum à la fréquence ©, La relation 


O3 T3 — À (8) 


est donc vérifiée. 

Supposons maintenant qu’à cette fréquence «, 
nous fassions une mesure en fonction de la tempé- 
rature. Avec les notations antérieures, €, a son 
maximum à la température T,, définie par l’équa- 
tion (7) ou sa transformée par l’équation (3) 
dérivée : 

aU 
RT2 


(1 — w2 72) — (1 + 2 +2). 


(9) 


Éliminant ©, entre (8) et (9) et mettant en évi- 
dence le rapport : 


E- 
T3 5 
on a successivement 
ie NE PEN 
es ll == =") 
RS [: Fe : E, Ta 
pis 


\ k TT.) [< Se 


avec la condition 


EN UREL  N aUIRT € 


ER Na 10 
T, T, 2U À (aUJKTi) + € (0 


si & — Cte on retrouve évidemment 7; = Ta. 


2.2. RELATION ENTRE TEMPÉRATURE ET FRÉ- 
QUENCE DU MAXIMUM DE €. — Sans vouloir dis- 
cuter ici des limites physiques de validité de la 
formule connue 


Vi rAen tre 


remarquons qu’elle est destinée à coordonner des 
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mesures faites en fonction de la fréquence’ et qu’elle 
doit donc s’écrire avec nos notations : 


Vs = À exp (— UJKT';). (11) 


Nous aurons sa nouvelle expression pour relier 
des résultats expérimentaux pris en fonction de la 
température, en y remplaçant T, par sa valeur. 

Cette expression en fonction de T, s’écrira sim- 
plement à partir de l’équation (10) (on peut au 
premier membre remplacer v, par v, sans incon- 
vénient). 


U 
Ie SE = le il — RT, 
DS SNaU AT ee 
log v,— - © 
Be Br eo NE UE + € 
ou 
_ QUIRTE) — & _yxr 9 
ne AVES Ua 
Si donc l’on reporte dans un diagramme 


(log v, 1/T) les couples expérimentaux (vs, Ti). 
relevés en fonction de la température, ils ne 
s’alignent plus sur une droite mais sont distribués 
sur une courbe à concavité positive ou négative 
selon que €, est négatif ou positif (fig. 3). 


log 


IDE Sh 


Il revient au même de considérer que la for- 
mule (11)-est valable aussi bien en variation de 
fréquence qu’en variation de température, soit en 
considérant que À est une fonction de la tempé- 
rature et posant : 


D MaUIREE) Ze; 
JS QUES Les 


où A, et U sont des 
constantes 


(13) 


soit en considérant que Ü est une fonction de la 
température et posant ; 


(aUJKT?) — & où Aet U,sont des (14) 
U = U; log VUE constantes. 
2.3. CALGUL DU MAXIMUM DE € — Il suffit 


d'éliminer wr et r’ entre les équations (7), (2, (3) 
et l’équation dérivée de (3). 

On obtient pour expression de cette valeur 
maximale : 


om = U J2U) Va? Ur — KT 6? (15) 
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dépendant de la température, de l’énergie d’acti- 
vation et de Es. On ne retrouve Em — Cte que 
dans le cas où &, — 0. 

On peut du reste prévoir à quelle température 
(et donc à quelle fréquence) l’ordonnée €,» sera 
maximale : il suffit de considérer l’expression (15) 
comme l'équation de la courbe enveloppe des maxi- 
mums de e, en fonction de la variable 7}. 

Il y aura donc «maximum maximorum » pour la 
température 7, telle que 


dEom [dT — 1) 


solution de l’équation 


aU? — k? T2 + Teï] = 0 (16) 
(nous avons supposé U — Cte). 
2.4. Discussion. — Nous avons montré que les 


couples (w,, T,) relatifs à un maximum de €, à la 
fréquence fixe ©, en fonction de la température 
étaient différents de ceux (w3, T,) relatifs à un 
maximum de e, à la température fixe 7°, en fonc- 
tion de la fréquence. En particulier, si l’on impose 
©, — &,, on a en général 7°,  T,, puisque 


(aU JRT) 
au JRT?) + &° 


té 
— E&s 


Suivant que s, soit positif ou négatif (proximité 

d’un changement de phase) l’écart (z — r. 

2 3 
change de sens, mais sa valeur absolue ne dépend 
que de la valeur absolue de ei et croît avec elle. 
Il tend également à s’accroître pour les grandes 
valeurs de T'et les faibles valeurs de U. 

Avec les données usuelles, la valeur absoluz du 
logarithme ne dépasse pas quelques unités, la diffé- 
RÉMETe 
maximal k4/2U. 

Par exemple pour T, — 200 0K et U = 0,1 eV, 
supposant le logarithme voisin de 1, on trouverait 
T; = 2220K ou T; = 182 °K suivant le signe 
pris pour €. 

Il est courant de calculer l’énergie U en reportant 


rence a donc pour ordre de grandeur 


dans un diagramme (ie V, 7) les couples v, et 7, 


relatifs à des maximums de e, pris en fonction de 
la température, et mesurant la pente de la « droite » 
passant au mieux par les points expérimentaux. 
Le résultat reste approximativement valable si la 
quantité 


(QaUJRTY) — € 
(aUJRTS) + es 


varie peu avec la température, mais il peut en 
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No 


être différemment et les points expérimentaux ne 
pas s’aligner sur une droite. 

Le calcul (formule 16) rend bien compte de 
l'existence de maximums de €; d'amplitude inégale 
passant par un maximum à une température 
donnée. Malheureusement, l’expression fait inter- 
venir la dérivée seconde de &,(7) qui n’est en 
général connue qu ’empiriquement. Il est donc pos- 
sible que les prévisions de cette formule soient 
purement qualitatives. 

En particulier, supposant une variation linéaire 
ou parabolique pour (T°) et faisant les hypothèses 
simplificatrices e,(0) = ex (cf. [B]) et e (0) = 0, 
on trouve deux valeurs 7, indépendantes des para- 
mètres de variation, soit : 


Tm = 8 200 U (variation linéaire) ou Tm = 4 739 U 


où Uest l'énergie d'activation exprimée en électron- 
volts. L'ordre de grandeur reste trop élevé. 


III. Vérifications expérimentales. — Les résul- 
tats de la théorie précédente seront d'autant plus 
faciles à vérifier que dans la gamme de température 
considérée, la variation de e,(T) sera plus grande. 
Il est également indispensable que le corps choisi 
présente des phénomènes de relaxation diélectrique 
dans la région explorée. 

Nous avons choisi les deux aluns d’aluminium 


No 


ammonium et de fer ammonium. Ces deux aluns 
présentent un changement de phase connu aux 
températures respectives de 55 0K et 85 0K. Au- 
dessus de ces températures, #; est négatif dans un 
intervalle assez grand et une relaxation dipolaire 
_y à été mise en évidence depuis longtemps par 
Guillien [6]. Les courbes expérimentales que nous 
avons reprises et confirmées présentent bien l’allure 
des courbes théoriques de la figure 1. Examinons 
séparément nos résultats pour les deux corps. 


3.1. ALUN D’ALUMINIUM-AMMONIUM. — Nos 
mesures, comme celles de Guillien, portent sur l’alun 
fondu, et sont faites en fonction de la température 
à des fréquences de 50 Hz à 1 MHz. A la précision 
des mesures, nous pouvons confondre €, avec 
1508) et Sa variation en fonction de T'est repré- 
sentée par l’expression 

960 
7 T+88 
entre 62 et 140 OK à moins de 5 % près (en utili- 
sant les valeurs numériques prises sur les courbes 
de Guillien). 

Nos vérifications ont porté sur les points 
suivants : 

1. On peut comparer les valeurs calculées et 
expérimentales de la différence & — €; à plusieurs 
températures. 

Le tableau ci-contre montre le bon accord entre 
les deux séries de résultats : 


Es 


Es Eym 
T°K MESURÉ CALCULÉ 
102 0,16 0,17 
95 0,15 0,16 
80 0,13 0,14 
65,5 0,10 0,12 


2. On peut.également rechercher à quelle fré- 
quence se produira un maximum de €; pour une 
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température donnée (équation 5). Prenant ainsi 
T, — 6595, le calcul fournit 


v = 5 670 Hz 


pour une valeur expérimentale de 5 700 Hz. 

La concordance est moins bonne aux tempé- 
ratures supérieures, ce qui peut provenir de la 
relation exponentielle entre T'et v. 

3. À partir de l’équation (15) nous pouvons 
calculer l'amplitude du maximum de €, connais- 
sant la température à laquelle il se produit. Le 
tableau ci-contre résume la comparaison entre les 
résultats du calcul et ceux de l’expérience : l’accord 
est moins bon aux températures élevées. 


ED , °2m ; 
MESURÉE CALCULE 
63 1,80 1,92 
70 1,63 1,67 
74 1,50 1,53 
78 1,27 1,41 
81,5 1,04 1,30 


4. Opérant comme pour €,, nous pouvons cal- 
culer directement à partir de l’équation (16) la 
fréquence pour laquelle e, passera par un maximum 
à une température déterminée. 

Nous comparons ces valeurs avec celles déduites 
des résultats expérimentaux par interpolation sur 
une droite (log v, 1/T), en utilisant les corrections 
du paragraphe 22. 


TK log v 
A 
48 T3 MESURÉ CALCULÉ 
65,5 (5,7 4,18 4,5 
75,5 76 5,42 9,7 
85 86 6,34 6,56 


5. Par contre, la température T,} déduite de 
l’équation (16) est d’environ 152 0K alors qu’elle se 
situe en réalité vers 65 0K. 


IBNES 


614 JOURNAL DE PHYSIQUE 


3.2. ALUN DE FER-AMMONIUM. — Les résultats 
expérimentaux sont résumés dans la figure 5 et 
confirment également des résultats antérieurs de 
Guiïllien [7], Freymann [4] et Powell [9]. Les 
courbes e,(T) et e,(T) ont une allure générale 
identique à celles de l’alun d’aluminium. Notons 
toutefois que pour les fréquences inférieures à 350 
kHzenviron, la température T, solution de l’équa- 
tion (5) reste inférieure à T,, d’où des courbes de 
permittivité réelle pratiquement confondues avec 
(T°) et passant par un point anguleux à T4. 

Aux fréquences supérieures, on retrouve le dépla- 
cement de Em. 

La courbe e,(T) peut se représenter dans l’inter- 
valle 84-150 0K par la relation : 


es = 26,5/VT — 73,5 


avec une erreur inférieure à 3 %, [3] ce qui nous a 
permis les vérifications suivantes : 

1. Dans le tableau ci-dessous nous indiquons 
pour 4 fréquences les températures 7, et les 
valeurs de e,, correspondantes déterminées par 
lPexpérience. 

Pour chaque température, nous avons calculé 
la valeur e;, à partir de (6) et la fréquence à 
laquelle devait se produire le maximum. 


TOoK y log Y Em 
KHZ MESURÉ  CALCULÉ MESURÉ  CALCULÉ 
91,2 500 5,70 5,74 5,60 5,71 
DRE AUDE 6 5.98 5,20 540 
99 2000 6,30 638 47% 490 
104,8 5 000 6,70 6,78 4,30 4,28 


2. Nous avons également comparé les valeurs 
calculées et expérimentales pour €,(7) à 1 MHz, 


à des températures inférieures à la température de 
raccordement avec £{T). 


En 
TK MESURÉ CALCULÉ 
85 3,38 3,20 
87 k,12 k,11 
89 k,54 k,74 
93 5,16 5,36 
95 5,18 5,38 
97 5,14 5,34 
99 4,96 DZ 


3. Aucune conclusion ne peut être tirée de la 
formule (16) puisque la température T,, inférieure 
à la température du changement de phase n’est 
pas observable expérimentalement. 


Conelusion. — En définitive, il apparaît donc 
qu’il est indispensable, lorsque des mesures de £e* 
sont faites en fonction de la température, de tenir 
compte dans les calculs ultérieurs de la variation 
de &, avec cette température. Les calculs nécessités 
sont plus longs et délicats que lorsque l’on opère 
en fonction de la fréquence, mais cette difficulté 
est largement compensée par les simplifications de 
l'expérience elle-même. En même temps nous 
pensons avoir donné une explication valable pour 
certaines « anomalies » souvent rencontrées dans 
les courbes expérimentales de constantes diélec- 
triques en fonction de la température : maximum 
de €,, amplitude variable selon la fréquence, des 
maximums des €,, courbure des lignes représen- 
tatives dans le diagramme log v, 1/T, désaccord 
subsistant entre différents auteurs. En particulier, 
il sera toujours indispensable, lors de mesures en 
fonction de la température, de chercher l’ordre de 
grandeur de la différence T, — T; (équation (10)) 
qui peut devenir notable. 


Manuscrit reçu le 20 avril 1960. 
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LETTRE À LA RÉDACTION 


ANALYSE DU COEFFICIENT DE CORRÉLATION 
DE LA FONCTION ALÉATOIRE X(#) — + 1. 


APPLICATION A L'ÉTUDE DE LOIS 
DE DÉSINTÉGRATIONS DE RADIOÉLÉMENTS. 


Par MM. G. Lanpaup et Cl. MagBoux, 
Laboratoire de Physique Corpusculaire de Caen. 


I. Introduction. — L'étude des lois de désintégra- 
tion de radioéléments peut s'effectuer de la manière 
suivante : À l’ensemble Æ des événements aléatoires 
6, 6541 … constitués par l'émission de particules, on 
associe la fonction aléatoire X(t) = + 1 [X(t) bas- 
cule alternativement de + 1 à — 1 à l’apparition de 
chaque événement 6;] [1]. 

L'étude du coefficient de corrélation 


CR) =E{f X@.X(+R)} 


en fonction du retard À permet de déduire la loi de 
désintégration du radioélément envisagé [2]. 


IT. Description de l’appareillage. — Quel que soit 


l'instant t choisi, le produit X(t). X({- k) ne prend 


que les valeurs + 1 ou — 1. Expérimentalement, il 


_ suffit donc de déterminer le signe de X({).X(t + h) à 


N instants quelconques f,, t, ... ti, ... tx et lorsque 
N est suffisamment grand C(h) est équivalent à : 
1e 
7 Dizs Au). Xi + h). (1) 


L'appareil permettant d'effectuer l'opération (1) 
est décrit par le bloc diagramme ci-dessous : 


Le détecteur D reçoit les particules émises par la 
source radioactive A. Les impulsions issues de D 
déclenchent le multivibrateur bistable B. 

A la sortie de B, on dispose de la fonction aléatoire 
X(t) = + 1 associée à l'émission de particules de la 
source radioactive À. 

L'ensemble, A, B, D constitue donc le générateur de 


la fonction X(t) = + 1. 7 
A l'instant t — t; le générateur G délivre une impul- 


sion brève, qui attaque le circuit de coïncidence C,. 
Cette même impulsion retardée par le circuit de 
retard R attaque en { — tj + h le circuit C, de concep- 
tion analogue au premier. 

D'autre part C, et C, sont simultanément attaqués 
par le signal X(t) — + 1 qu'il faut analyser. 

Le circuit de coïncidence C,; transmet l’impulsion 
brève issue du générateur lorsque X(t;) — + 1, par 
contre C. est bloqué si X(t) — —1. 

Le circuit C,, identique à C;, transmet l’impulsion 
d'attaque si X(t4; + h) — +1. Enfin le circuit de 
coïncidence C; délivre une impulsion lorsque les deux 
conditions suivantes sont vérifiées. 


19 ÆX{tiÿ) = +1, 
20 X(tj + h) = +1. 


Cette impulsion finale est enregistrée sur l’échelle 
de comptage S:. 

Pendant le temps T, le générateur G délivre N impul- 
sions, qui sont enregistrées par l’échelle de comptage 
S;. S, indique alors n événements qui vérifient les 
conditions (2). LE 

Lorsque N et T sont suffisamment grands, le 
d'événements vérifiant : 


X(tiÿ) = —1 
At + h) =—1, 


(2) 


ombre 


tend vers n si E { X(t)} = 0. Donc, parmi N expériences, 
il y en à 2n pour lesquelles 


At). Ati + k) = 


Le coefficient de corrélation 


44. 


CE) 7% DIET Alu ).Xlt + ho), 


tend vers 
Ch) —1/N.[2n 


(N — 9n)] = 4n/N —1. 

III. Domaine d’utilisation et essais du corrélateur- 
— Le circuit de retard R(h) permet l’analyse continue 
de C(h) pour des retards variant de 3.105 $s à 
DAV, 

L'activité de la source A est déterminée de manière 
à ce que le nombre moyen pÇ d'impulsions détectées 
par seconde soit compris entre 10% et 2.10%. Dans ces 
conditions les circuits de coïncidence C, et C, sélec- 
tionnent les événements ‘X({) — + 1 à mieux de 
2 0/60- 

FF d'évaluer la précision des mesures, l’essai sui- 
vant a été effectué : On a substitué au détecteur D un 
générateur d’impulsions périodiques de fréquence de 
répétition po — 5 000. Le coefficient de corrélation 


616 


C(h) correspondant à cette distribution périodique 
d’impulsions est égal à : 


CR) = (1) (2n + 1 —2p6 7) 


ET D EURE E 115 


Cette courbe théorique est représentée sur le gra- 
phique I. L'écart avec les points expérimentaux est 
inférieur à 3 2/60. 


GRAPHIQUE I. — Expérience effectuée avec : 
00 === 750, 
Po = 1 500, 
0610007 


IV. Résultats. — Nous avons finalement analysé 
une distribution poissonnienne stationnaire d’impul- 
sions due à l’émission de particules d’une source de 
Coso- Le calcul de C(h) donne C(h) = e—28h [3] (06 
étant le paramètre de la loi de Poisson). 

Les courbes expérimentales C(A) (graphique IT) ont 
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a Expérience pour @, = 2500 
» n= 4000 
»n = 8000 


GRAPHIQUE Il. 


été obtenues pour diverses valeurs de p, (2.10, 
41098108) 

L'écart entre les points A et expérimentaux 
est inférieur à 5 2}, lorsque p5 h < 1. 

Par contre, la précision est de l’ordre de 1 % lorsque 
LORS 

Les résultats obtenus permettent de penser que la 
méthode décrite peut être utilisée pour déterminer des 
lois de distribution non poissonnienne ; en particulier, 
pour des corps radioactifs à filiation, C(k) est essen- 
tiellement fonction de la période de l’état métastable. 


Lettre reçue le 10 mai 1960. 
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